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Formulas for Inverse Osculatory Interpolation 


in the Complex Plane 
Herbert E. Salzer ' 


Improved formulas for inverse osculatory interpolation in the complex plane are ob- 
tained by inversion of Hermite’s formula and the use of appropriate grid point configurations. 
They cover the cases for n=2(1)7, where n is the number of points required in direct oscu- 
latory interpolation. The formulas provide an improved means for inverse interpolation 
in the complex plane where the first derivative is either tabulated alongside the function 
or is easily obtained. 


Formulas for n-point inverse osculatory interpolation for finding s=2)+ ph from f(x) in 
cerms of f,=f(a,) and f,=f' (a), where the ,=2)+kh, the k ranging from —[(n—1)/2] to 
(n/2], are equally spaced at intervals of h, have been obtained previously by the author [1] ? 
through the inversion of the Hermite osculatory interpolation formula [2,3]. Those formulas 
are still applicable to analytic functions when z is replaced by a complex variable z=2,+ Ph, 
and 2,=2)+kh are equally spaced upon any line in the z-plane at interval h, which may be 
any suitably chosen complex number. But for many analytic functions that are given with 
their first derivatives over a Cartesian grid of length A in the complex plane (so that now h 
is a positive number), the accuracy of inverse osculatory interpolation is considerably greater 
when the points 2, are chosen as grid points that lie much closer to each other, no longer in a 
straight line bevond the 2-point case. 

Improved formulas for inverse osculatory interpolation, which utilize the Cartesian grid, 
will be particularly useful in connection with such tables for complex arguments as (1) log ['(2), 
together with its derivative the psi function, (2) Bessel functions of the first or second kind 
giving Jo(z) and J,(z2) Ji(z), Yo(z) and Y\(z)=—Y%(z), or linear combinations of them, 


(3) probability integral e~“duwith its integrand e~* [4,8], (4) miscellaneous tables of integrals 


of the more elementary functions where the first derivative or the integrand, although not 


tabulated, is easy to calculate, viz., the function (e~“/u)du [5], and (5) tables of solutions of 
important linear differential equations, together with the first derivative [6]. In all such tables 
and in many others, the user will find the inverse interpolation formulas that are given below 
to be particularly convenient, especially in those cases where grid length h is too large for 
sufficiently accurate complex inverse interpolation, using the formulas in terms of only the 


functional values f; {7}. 
In accordance with the choice of grid points that has been employed consistently in earlier 


papers of the author on complex interpolation, for n ranging from 2 to 7, the n-point formulas 
in this present article are based upon the values of f(z.) and f’(z,) for z, in the following 
configurations: 

Two-point Three-point Four-point Five-point 


Six-point Seven-point 
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These formulas enable one to find z=2)+Ph by expressing P in terms of /(29+ Ph) =f. 
tz) =fe, and f’(z,) =fi, for 2,=29+kh, where k ranges over the points of the selected n-point 
configuration, the choice of n depending upon the number of points required in direct osculatory 
interpolation. Although the direct interpolation formula for n=6 and n=7 is of the lith- 
and 13th-degree accuracy, respectively, (recalling that the n-point direct osculatory interpola- 
tion formula is of (2n—1)th-degree accuracy), the inversion formula for ? given below does 
not go beyond the 10th-degree terms. In fact, most practical problems will require only the 
first few terms, and it is very rarely that one would need this inversion series up to the 10th- 
degree terms. Thus in this paper we do not use the coefficient of P™ in the 6-point Hermite 
formula, nor the coefficient of P", P'’, and P® in the 7-point Hermite formula. 

For every n we define r= (f—fo)/hfo, and corresponding to each of the above n-point con- 


figurations, where n=2(1)7, quantities s, ¢, u,v, w, 7, y, 2, and Z are defined as follows: 


Two-point 
s 3{ fo—hi} —h{2fotsi})/hfe, 
t=(2{ fo—fis +hi fotty)/Afe 
u ? u z y 2 0 

Three-point 


s= }(8ifp+ (1—42)f,—(1+-42) fh, +h (4—47) f,+7f1—fi})/hfo, 
t= 4(— (12+ 122) f+ (7 +-52)f, + (54-70) f, +h{ —8ify— (2 +a fit+ (2—a)f;} )/hfo, 
u=4(14f.—(7—27)f,— (7 +27) f, +h (44+-42)f6 + (2—1)f 1, — A — 22) f3} )/hfo, 


v= 4(—(4—41)fp+- (1—32)4,4+-(8—af,+-hf —2f04+0f;—7f;})/hfo, 


wx s=y= s=72=0 
Four-point 
s=4(197f,— (8 +-62)/, + (8S—62)f,—Tifi., +h (6—67)f,+2f,—27f;— (U1 —a)fis:})/Afo, 
t= }3(—(374+-377)f.+ (382 —127)f,— (12—322)f, 4+ (17 +- 170) fini +h! —17if,— (4—60f14+ (44+ 60)f; 
5if 1, ) hf, 
u=4(70f.—(15—502)f, — (15 +-507)f,—40f,, , +A! (144+ 147) f,— (44+-112)f,-—(114+-40)f; 
+-(6+-67)f 14.) /Afo, 
v= }(— (388—381)fp— (25+:371)f, + (37 +257) f, + (26—262)f,, +h! —14f)4+ (10430), 
(10—32) fi—S8firi) )/hfo, 
w= }(—23if,+ (21+4-22)f, — (21 —27)f,4- 197; hi (4—47)f,—(4—37) f, —(8—40)/f; 
(3—32)fisi} /hfo, 
r= $((38+-32)fp— (8—32)f, + (8—32)/ 3+31)f hi if,—ih,;—it +f) hf, 
y=n = 3=0. 
Five-point 
s= (zi, (1875 —62507)f, — (4000 +- 60007) f, + (108 + 2547) f. + (2272 + 4967)f 


+ (3500 — 10002)f, , (200— 1507) f, + 200f, — (3 +-42)f, + (82 — 247); 


50 
(50-4 502) f; hf, 


234 


tl 


1 (1525 4+-39507)f, + (5600 + 12007) /, — (223 +-1147)f,— (1952 —464z)f, 


ivo 


{ 


5 { (65 —23072)f,— (240—2407)f, + (9 +27) f,— (32 —641)f; 


(1900 —24007)f 


 (120+-202)fi4:})/hfo, 


(salvo | (59750 + 371252) f,— (96000 — 530002) f, + (5284 — 13132) f2 + (29216 —263122)f, 

L (1750—625007)/, + hs { (625+ 16257) f, + (100—34002)f; — (89 — 487) f, — (24 +7322) f; 
(1125—5752)fis4})/hfo, 

(sang { — (72625 — 52507) f. + (48000 — 1180007) f, — (4103 —57467)f.— (14272 —435041)f, 

- (43000 +-635002)/, sh { — (2800 + 18257) f, + (5200 + 62007)f; + (96 —2472)f; 

- (896 +-15287)f; + (1450 —28007)f,,:})/hfo, 

(yh e{ (3850 —30757) fy + (2450+ 92007) f, — (42 +5817) f.— (658 + 33442)f,— (5600 4+- 22002) f, 41} 


n{ (215—52)f,— (600 4+- 1102) f, + (7 +212) f,— (1144-722) f;, + (304-2502) f;, «})/Afo, 


(obs 3500 — 101257) f,— (17250 +- 1475072)f, + (1022 + 11217)f, + (4978 +-6254z)f, 

- (14750—27507)f, cpg —(875—2502)f, + (1300 —6507)f,— (49 + 322)f, + (316 +-382)f; 
(400 +-5002) fi.) )/Afo, 

= (1000 +- 53757) fo + (12500 + 20007) f, — (898 + 1897) f.— (3852 +- 19367); 
(6750—55002)/, rpg (50—1007)f,— (200 —3507)f, + (18—2)f,— (82 — 242); 
(1754-75) fis.) /Afo, 

(ralpo ( (375 + 5007) fp — (1500 —5002)f, + (117 — 442) f2 + (508 +-442)f, + (500— 10002) f, 4. ¢} 
gg { 25if,— 100if, — (4 —32)f.+ (16 — 127) f;—50f;,:})/hfo, 

0. 

Six-point 

(yop, 165007/9— (14336 + 11522)f, + (289— 9872) fo + (14336 — 11527) f,— 140007, ; 
(289 +-987)f2 -% (400—4007)f) + (192 — 2567)f, — (4—32)f34+- (256 —1927)f; 
(200—2007)f; (3—47)f3,})/hfo, 

(adie (86500 + 865007) fo + (91392 — 1026567) f, — (1219 —29837)f/.— (102656 — 913922)f;, 

-- (96000 +-960007)f (2983 —12197)f2;} +-525{ —31007/, + (896 + 30727)f;+ (4—532)f; 
(896 —30727)f; —2800%f; (4+-5372)f3;})/hfo, 


(ghy{ 95150fo+ (18688 + 1244167)f, — (1463 + 29347)f.+ (18688 -1244167)f,—129600/,,. ; 


(1463 —29347)f, *.{ (1500+ 15007)/5— (2816 +-11527)f;+ (38-4 397)f,— (1152 +-2816:2)f; 
- (1920 +-1920/)f; 39 +-387)f5,})/hfo, 
(rasa { — (169500 — 1695007) fo— (301824 + 1927687)f, +- (9957 4+- 26017) f.+ (192768 +-3018247)f, 
-++- (271200 —2712007)f, (2601 +99577)fo;} +425{ —4970f5+ (7808 —42247)f; — (172—%)f; 
+- (7808 +42247)f' —8160f; (172+-92)f3,})/hfo, 
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10477507/y + (1662976 — 4663687) /, — (44774 — 3074387) f.— (1662976 + 4665687) / 


w= (aden 
1919000z/, (44774 4-30743))f, iy (14650 — 146507) f, — (8672 — 424967)/ 

(564 —7237)f3— (42496 —86727)/ 29300 —293007)f; 723 —5647)f5,))/hf 

r= (yaebwa) (919000 — 9190007) f, 1038592 — 208025672) 7, + (8369 — 790337) f. 

t- (2080256 — 10385927)/ (1890000 + 18900007) (79033 — 83697)f, reo0. 248500f, 
(36096 + 462727)f) + (271 + 18287)f 36096 — 462727) f; + 586007 | 
(271 —13287)f5,!)/hfo, 

y= (ats 28450/)— (19568 + 5O09767)/ 1393 + 14497)/,— (19568 — 509767)/, + 64800/ 
(1393 — 14497)fo;} +-y25{ — (750+ 7501) f, + (27524 1447)f 61+337)/3+ (144 +2752/)f 
(2040+ 20407)f; (33 4-6172)f5,')/hf 

(yalgg | (12000 — 120007) fy + (82224 + 127687)f, — (1357 — 997) f, — (12768 + 322247) 
30000 — 300007) f (99—1357,)f. 25 / 310f5— (608 — 6247)f; + (22—97)/ 
(608 +-6247)/; + 960/; 22 +-97)f5,})/hfS, 

2 cpg | 157500, $0688 — 147847)/ $12 — 11597) f, + (30688 + 147847)/,—430007/ 
(2124+-11597)f =m 200 — 2007) 64 -+-S487)f' 7 —247)f5+ (848 + 647)/ 
700—7007)f; (24—7r)/ hf 

Seven-point 

‘ 1060 (1960—21780r)/ 48400 — 112007)/ L119+ 15427)/, 15744 + 145927)/ 
28000 — 460007)/ 2432 —647)f 167 —2947)f, —_ {80 —4407)/ 

SO07f, +- (24 +-77)f,4+- (384 4-1127)/ 1400+ 2007)f, 32+-162)f; 3-+40)f- hf 

t= (sek; (1133164 1470127)/ 193280 —4070407)/ 21313 + 16597) f, — (242688 + 307847)f; 
120000 + 6000007) / 20736 — 209287) 675+-41057)/ soo. (344—40927)/ 
7040 +-60807)/ 268—1517)f 1448 — 15367)f; 13600 — 92007) $16—1287)/ 
(60 +-152)f,,) hf, 

U= (qo'eu | (825894 + 1402087 )/ SYS 176 + 228405327)/ 94875 —832507)f, 

t- (1152000 —6720007)/ 2T7HO0000 — 19920007) (2688 728007) f, 
(18507 4+-162907)f, og { (10044 + 129087) f,— (51392 — 66567)f, + (470— 17652) f- 
- (14880— 217607)/; — (21600 — 904007) f; 1248 —22407)f,- (333 —1587)f3,})/Afo, 
v 5006 (2094204 — 13613727)f, 7780992 + 32285447)f, + (24441 —46638772)/ 
(1650816 — 42565127)/ 11844000 — 14600007 )/ (433920 + 4663047)f, 
+- (91491 4-335 12)f, se (44834 +-74887)f) + (103744 — 1473927)f, + (3484 4+-58377)/ 
(4704 + 918727)f,— (160800 + 2896007) / 2400 +-91847)f,. ,4+- (588 —12457)f5,) /Afi, 
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= (662260 — 25734307) f, + (9236480— 39006402 hi + ( 104714 1376775 fs 
(1662464 51515527)f (9358000 1107 10007)f, (8107524 164961) fo, 
2112807)/, — (8494 4+-20177)f; 


93742 —704497)f, seo. (88730—251907) fo + (43840- 

71904 +-836727)f, 387500 + 1301007)f; (10352 + 63767) f; (423 +-17247)f,})/hfo, 
(paboo (1955376 + 36844327) f, — (6505920 — 159945607), — (1095693 + 279997) fo 

7805568 + 50226247)f,— (1074000 + 254340007) f, (1139136 — 9698887) f,, 


53805 — 2095050 )f revo, — (17684 — 679627) f,— (3153604 1971207)f; 
13623 —S86867)/ 186528 +371047)f; — (667600 — 2930007)f, (22016 —46887)f, 
26154-20102). hf, 

. 1149156 + 4146427)f, — (2086704 + 5056 1287)/, + (264015 — 2577307) 

2906640 — 3424807)/; + (5796000 + 58230007) /;, + (52992 5037607) fs 

29493 + 667807) f, “= (9306 + 177927)f, + (114168 —278247)f, — (1095 — 53657) f, 
53460 — 273207)f; + (98100 —2179002)f; , ;— (4332 — 63202) f)4 , + (1122 — 127) fs} )/hfo, 


985056 — 4008087) f, + (4675648 + 20239367) f, + (221 4+-3634037)f, 


1877904 — 19585287)/, — (7782000 + 3940007)/f, (306480 + 4002567)f/2,, 


932487) f\— (30714-44787) f; 


63309 + 342597) f, revo. (16076 + 59827) f,— (56336 
23024 — 505687)f; + (69600 + 2194007)f; (1200 +-75767)f; (877 —6907)f3, })/Afy, 

== 74180— 1509907) f, — (802080 — 2974407)f, — (45912 + 448917) /.— (85088 —4520167)f 
1034000 — 8390007) 89296 +- 146087)fo., + (12556— 19477) fos) + she (2490—9707)f, 
1240 + 176807)f, + (952+. 1617) f, + (2232 + 89767); — (343004. 205007) f,, , + (1096+ 8487)f,; 
D1 1927 I. Af 

For every n, P is given by the following formula: 

rs +7°(2s / y' 5s 5st — uw) +7°(148*—2187t4+- 3t° + 68su—r r $2” +-848°t —28st- 


7 ’ ‘(1328°—330s8t+ 180877 1L2Os°u 127 T2 stu — 36s8*0 + 47 + Str 


IN Ttu rey, u } 2 
Rs —a rs 129% L2878°t —9908°¢ $95s8'u $95s-tu 165s/ L658°7 t5t°u $5s8u- 
Odes, $5s Ww iy) Of Qeyy y + 7 ( ] L308" 5OOS8t + SOOSS8 4 2OOYs u ] 1308-7 


PROOS Tu 1158*¢ jdt HOOst-u 3308-u- H60s87tr | 2POs?w— 5dtu jotr 110sur 
avp LOww LOta LOsy : ,* N62" 1944887 — 2402485? — SOOS88°u 


L1Ostw s+ >! 


LOOLOSs°? L5OLSs*tu —300388°7 LOOLst?— 60068777 u — 2OO2S*u $O04s°t7 LOOLs*w 
ISO U S5OSstu S5S8s-ti + S5Sst-7 + S5ASs-ur —2S68° 2 aa} L32tur — 660 w— 6680" 


132suw— 132sta—6682y+ Llew+ Llua-+ 1ity+ 1182 
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1957 Research Paper 


An X-ray Study of Textural Stresses in Two-Phase Alloys 
C. J. Newton and H. C. Vacher 


Internal or residual stresses of the textural stress class are discussed, with special atten- 

tion to that type resulting from a difference in thermal expansion of two phases present in 

in alloy X-ray stress Measurements were made on various brasses and steels, and the 

esults were found to be in reasonable agreement with the quantitative theory of F. Laszlé. 
1. Introduction were assumed to have surfaces free of normal stress 
and to pack together to form the body of material, 
The term “textural stress’? will be used in this although strictly speaking only prismatic or poly- 
dy in the sense adopted by E. Orowan in a sym- | hedral shapes could do so. The formulas for the 
nosium of the Institute of Metals in 1948 [1 His | approximate values of the principal stresses in these 


most inclusive term, internal Stress, is defined to he 


stress existing in a body upon which no externas 


force Is acting Some writers use the term ‘residual 
stress” in this wide sense This field is divided into 
two classes One class consists of body stresses 
which acise when an external factor has affected 
different macroscopi1e parts of a bods differently, 
} even though the material of the body may have been 
quite homogeneous The other class consists of 
textural stresses, which result from inhomogeneities 


activating external 

influence may have itself been homogeneous. Under 
| ordinary circumstances, where the grain size is small 
compared to the MAacroscopt1e dimensions of the body, 
equivalent the widely 
ised terms, Macrostress and microstress, respectively. 
The former implies a system constant over 
many contiguous grains in a body, balanced by 
opposing stresses In other grains at some distance; 
the latter implies a stress system that varies from 


of the material, even though the 


the classes are closely to 


se 


st ress 


grain to grain or within a grain, and is balanced over 
short ranges. 
The most detailed) and quantits ative theoretical 


analysis of textural stresses arising from the presence 
of two phases in an alloy is to be found in a series 
of four articles between 1943 and 1945 by F. Laszlé 


2|. The four articles were reviewed in a nonmath- 
ematical manner by F. R. N. Nabarro in 1948 [3]. 
Laszlé discussed textural stresses arising from 


various sources, such as anisotropy of the grains and 
differences in thermal or mechanical properties, and 
the effect of external stress on the textural stress. 
Much numerical detail was presented with regard 
to iron alloys. The greater part of the work dealt 
with an involved mathematical analysis of textural 
by the difference of coefficients of 
thermal expansion of two phases present in various 


stresses caused 


proportions, and with various grain shapes. The 
solution was only approximate, of course, because 
of the simplifications that were necessarily made, 


In particular, 
were assumed: 


three types of idealized grain shapes 
Lamellar, with alternating lavers of 
the two constituents, and compound spheres and 
compound cylinders, with and of 
the two different phases. These compound blocks 


TT 


Figures in brackets 


covermgs cores 


ndicate the literature references at the end of this paper 


units presented in part I were used in this study. 
In Lasz'd’s article, part II, the ideas were extended 
to other special cases, such as capped evlinders, and 
the physical rigor was somewhat enhanced; but the 
results there were presented in the form of systems 
of simultaneous equations of high order, and the 
value the much more difficult solutions of these 
equations was questionable in the light of the experi- 


o! 


mental margin of error anticipated in the study at 
hand, 
Laszlé6 suggested four technological effects in 


which such textural stresses might be of importance. 
Yield, the formation of cracks, fatigue, and corrosion, 
the Jast three of which are most likely to be signifi- 
cantly influenced by such stresses. Although the 
mathematical erudition and the extensive calcula- 
tions from theory displaved in these articles were 
very impressive, no attempt at verification by X-ray 
measurements was made. On the contrary, it was 
specifically stated by Laszl6, part I, that these 
stresses could not be revealed by X-ray diffraction. 
In spite of this opinion, one might observe some 
stress with X-rays for the following reasons: X-ray 
diffraction is a highly selective process, one aspect 
of which is the sorting out of diffractions from crystals 
of the different phases at different Bragg angles. 
An averaging of the triaxial stress state over many 
grains of one phase within the material would result 
in a pseudohydrostatic stress. The relaxation of the 
component normal to a free surface, however, would 
result in a biaxial state of stress, the average value 
of which would be measurable by the two-exposure 
X-ray method described later. Even if the sensi- 
tivity of this method, which depends on the observa- 
tion of very small shifts in the diffraction line, were 
not sufficiently high for application, thermally in- 
duced textural stresses might be indicated by line 
broadening. In fact, broadening from this source 
was observed in a qualitative degree with duplex 
aluminum-silicon alloys by J. P. Nielsen and W. R. 
Hibbard in 1950 [4)}. 

If the experimental technique is skillful enough, 
however, actual net shifts of the line can be measured; 
such methods have been used by several investi- 
gators to study textural stresses due to differences in 
mechanical properties of two phases in plastically 
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deformed specimens. One of the earliest of these 
was D. V. Wilson [5], who in 1950 attempted meas- 
urements on the ferrite and cementite constituents 
in 1.26 percent carbon steel compressed 90 percent. 
Although his sensitivity of measurement of the 
cementite lines was very poor, he found a definite 
indication of opposing stresses in the two phases. 
In 1955 L. Reimer [6] published the results of work 
with iron and steel specimens deformed plastically by 
bending. He showed the influence of textural con- 
straints between the ferrite matrix and imbedded 
cementite. There was some indication that the 
residual stresses were proportional to carbon content 
by volume. In the same year, V. Hauk [7] made a 
study of plastically extended specimens of some two- 


phase alloys, not only steel but also a copper- 
aluminum alloy. He observed that not only an 
orientation dependence, but also the difference in 


elastic limits of the matrix and the intermetallic 
phase, contributed to the resulting stress system 

In view of the success in detecting by X-ray clif- 
fraction textural arising, above, from 
differences in mechanical constants in plastically 
deformed two-phase alloy specimens, it appeared 
worthwhile to investigate the possibility of observing, 
in annealed specimens, thermally induced textural 
arising from differences of coefficients of 
thermal expansion. It has been previously observed 
at the X-ray metallographic laboratory of the 
Bureau that textural are not completely 
relaxed by sectioning procedures [8]. For this reason 
it was hoped that the usual method of X-ray stress 
analysis on a free surface would permit detection of 
the strain resulting from a biaxial stress condition in 
the surface, which was up when the quasi- 
hydrostatic system of textural stresses in the two 
phases was sectioned. This is to say that, rather 
than the qualitative and somewhat ambiguous line 
broadening observed by some, the effect to be looked 
for was a small but dedinite line shift, which would 
indicate a nonzero net change in the lattice parameter 
of a given phase. The materials chosen were a two- 
phase brass and a 1.05 percent carbon steel. 


stresses as 


stresses 


stresses 


sel 


2. Stresses in Brasses 


Two types of brass were studied; one was a two- 
phase alpha beta brass, and the other was a single- 
phase alpha By adding a sufficient amount 
of zine to molten cartridge brass, a two-phase alloy 
was produced, containing 58 percent of copper and 


42 percent of zine. This alloy consisted of approxi- 


brass. 


mately equal volumes of the alpha and the beta 
phases. 
\fter a procedure of heating and rolling, small 


pieces of the proper size for mounting were cut from 
this material, and a heat treatment to give an appar- 


ent optimum combination of small grain size and 
minimum of preferred orientation of grains was 
determined. This was found to be 's hr at 900° F 
in air, followed by cooling in the furnace tube re 


moved from the heating coils 
The single-phase alpha material, cartridge 
brass, contained 70 percent of copper and 30 percent 


brass 


of zine. The absence of the beta phase was COD. 
firmed by metallographic examination and Xo 
diffraction. An optimum annealing condition ¢ 
this material, 800° F for }: hr, was similarly 
termined. 

The technique employed for finishing the Surfagy 
on the brass specimens consisted y many cycles a 
alternate mechanical polishing and elee tropolishing 
ending with the latter. The progress of the pro. 
cedure was followed by the growing sharpness of thy 
diffraction spots on stationary X-ray patterns, 

The diffraction patterns were made with X-ray 
from a cobalt target in the usual manner of the high. 
angle, or back-reflection, powder technique. 7, 
correct for changes of dimensions of the film and 
prov ide a precise value for the specimen- to-film 
distance, a calibration ring, caused by diffractig, 
from aluminum powder, was put on the film. Thy 
Bragg angle for alpha brass, using the (400) diffrge. 
tion, was about 75°, and for beta brass, using the! 
(310) diffraction, it was about 73 Two patterns 
were made with each specimen, one with the X- “ray 
beam normal to the surface of ithe specimen, and on 
with the beam making an angle of 30° with th 
normal. Only the second side emulsion was de. 
veloped on each film, which resulted in sharper line 
and greater contrast than is ordinarily obtained 

Measurements of the patterns were made with 4 
modified traveling microscope calibrated to the ord 
of 0.001 mm by the Shops Division of the National | 
Bureau of Standards. The Aa and Kaz lines wer 
in all cases, resolved, and in most cases readings were | 
made on each. Two complete sets of data wer 
taken independently by two different observers, wit 
results in good agreement. Simple precision analysis 
of the data indicated that the probable error in th 
average of each setting was than mm 
which, in a typical case, would correspond to 0.00008 
A in d-spacing. 

From the radi of the diffraction rings and th 
known Bragg angle of diffraction for the calibratior 
powder, it was possible to calculate the cosine of the 
complement of the Bragg angle of diffraction from 
the specimen. From this and the known value of the 
X-ray wavelength, the lattice spacing was calculated 

v Bragg’s Law. The two values of lattice spacing 
and the polyerystalline elastic constants enabled on 
then to calculate the by means of the well- 
known formula: 


less 


stress 


Ie l dt d 
l-++-p sin*y d 


where Poisson’s ratio, is 


is Young’s modulus, v is 


the angle between the surface normal and the dif- 
fraction plane normal, and the d’s are the latte 
spacings [9]. 

In connection with this formula, some investigators 


mechanical values ot 
Poisson’s ratio were satis- 
the values for beta 


it Was decided to 


have found that the ordinary 
Young’s modulus and 
factory [10]. However since 
brass were not well established, 
calculate these parameters for the alpha and beta 
phases from the fundamental single crystal elasti 
constants |1 1 by one of the polvervstalline agere- 
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vate-averaging procedures. The three principal 
- aedures were discussed in 1955 by H. H. Stadel- 
maier [12]. One formula, by A. Reuss, is based upon 
‘he simplifying assumption of equal stress in all 
one by W. Voigt assumes equal strains in all 
vrains. Of these two, the former has perhaps been 
ysed more, but G. B. Greenough [13] found that an 
average of the two results was preferable in his 
X-ray stress measurements. The averaging pro- 
cedure of A. Huber and E. Schmid [14] makes no 
assumption except that all grains behave as if they 
have equal elastic moduli, This procedure, though 
1 best agreement with observed moduli, is seldom 
ysed, perhaps because the answer does not proceed 
from a simple formula but from the numerical eval- 
yation of the leading terms in an infinite series, 
representing the numerical solution of an integral. 
Nevertheless, it was this procedure that was used to 
obtain the values employed in these calculations. 
These values were, for alpha brass, #=14.5 10° 
lbin2 and y=0.43; and for beta brass, E=9.95> 10° 
lb in? and y=0.29 
The observed relative lattice changes, 7 
d,, and the calculated values of the 
5 are presented in table 1 A positive sign indicates 


proced 


grains ; 


(dy 


st resses, 


TABLE | Stresse n brass 


tension and a negative sign indicates compression. 
This table also shows the probable error in the stress 
measurement and the theoretical values of the stress 
calculated for the two-phase brass with the formulas 
derived by Liszlé 

The value shown for probable error was based upon 
uncertainty in the average of the 
settings of the traveling microscope on the lines of 
the pattern of I< 10-2 mm, and an estimated 
incertainty of the temperature of +3° C during each 
of the two exposures. These estimations were prob- 
ably somewhat greater than was actually the case. 
They contributed the same order of magnitude of 


two sourees, an 


error, and, when they were combined, the total 
estimated probable error in stress was 3.210 
lb in.” This is in the neighborhood of the precision 


ordinarily claimed in X-ray stress measurements. 


because the Bragg angles of diffraction for both 
brasses were less than 76°. 

The theoretical stress values were computed with 
the aid of the formulas in the first paper of Laszlé 
previously referred to. It is to be expected that the 
stress revealed by the measurement of lattice spac- 
ings in the randomly distributed crystals of the 
aggregate would be proportional to the cubic dilata- 
tion of the grains. If & is the modulus of compres- 
sion, or bulk modulus, and p is the effective hydro- 
static pressure, the cubic dilatation is 


A pk. 
But p is related to 0, the first stress invariant, by 


p=-—9 a 
where 
0 o; YT Go-7 03, 


where o;, ¢, and o; are the principal stresses. There- 
fore 


A ( 
(1/k)-(0, + 02+ ¢3)/3. 


That is to say, the dilatation is proportional to the 
average of the three principal stresses. Therefore, 
it was this average stress, as calculated from Laszl6’s 
formulas and multiplied by a surface relaxation 
factor, that was compared with the observed value. 
The relaxation factor referred to is 1 minus Poisson’s 
ratio, according to Laszlé, equaling, in most cases, 
approximately 0.7, which corrects for the increase 
in stress components in the surface when the normal 
component is reduced to zero. 

The coefficients of linear thermal expansion used 
in the calculations were, for alpha brass, 22.4> 
10-°/°C, and for beta brass, 24.7 107°/°C. These 
values were calculated from data published by 
P. D. Merica and L. W. Schad [15]. The tempera- 
ture change was 455°C. It can be seen in the table 
that the results of the calculations of stress in the 
two-phase alloy were well within the experimental 
margin of error of the observed X-ray stress values, 
with the beta brass, having the larger coefficient, in 
tension, and the alpha brass, having the smaller 
coefficient, in compression. 

The value of stress measured for alpha brass was 
somewhat closer to zero than the estimated probable 
error, but the authors felt that, since it was carefully 
checked by both remeasurement and reannealing, 
this stress was at least qualitatively real. Although 
no definite second phase with a specifically differing 
thermal expansion coefficient was present, some 
thermal stresses could theoretically be caused by 
small amounts of disordered material of different 
thermal properties present at or near grain bound- 
It is possible, also, that some small degree 
might still have 


aries. 
of surface preparation distortion 


Although the actual error was probably not quite | been present, thus contributing to the nonzero 
this large. one could not expect oTeat precision result 
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3. Stresses in Steels 


The material from which the steel specimens were 
prepared contained 1.05 percent carbon, with minor 
constituents totaling 0.74 percent. Three specimens 
were prepared with different heat treatments in such 
a manner that one showed a structure of fine lamellar 
pearlite, one had a spheroidized structure, and one 
was decarburized. The decarburized specimen, con- 
taining a minute amount of cementite to 
saturation of the ferrite with carbon, was included 
to provide a specimen of an essentially single-phase 
alloy having the same composition as the ferrite in 
the pearlitic and spheroidized steels. Such a speci- 
men should contain no textural stresses ascribable 
to the presence of a second phase, and might proy ide 
a base line for the measurements on the pearlitic 
and spheroidized steels. 

The decarburized specimen was prepared by 
cutting \¢-in.-thick disks from the stock material 
and decarburizing them with wet hydrogen at 1,650 
F. Several attempts were made before a homog- 
enized specimen containing a trace of pearlite was 
obtained. The heat treatments to obtain the pearl- 
itic and spheroidized conditions were conventional. 
The decarburized portions were removed and the 
whole specimens were given a prolonged electro- 
polishing treatment in order to diminish the possi- 
bility of body resulting from a 
concentration gradient. The specimens then were 
mounted in Bakelite and polished by repeated eveles 
of light mechanical and electropolishing treatments 
until the sharpness of the X-ray diffraction indicated 
that the surfaces were essentially free of cold work. 

The same two-exposure method of obtaining X-ray 
data for surface stress calculations was followed with 
each of the three types of steel as had been followed 
with the brass specimens. From the steels, diffrac- 
tion lines were obtained only from the ferrite and, 
during the same exposure, from a calibrating powder 
spread thinly on the surface of the steel. Diffraction 
from the cementite is not observed in the high angle 
region from steel with this carbon content. Three 
pairs of patterns were made from each of the three 
specimens with three different radiations and slightly 
varying conditions. In one radiation 
was used with tungsten as the calibration powder, 
iron radiation was used with tungsten calibration 
also, and chromium radiation was used with silver 
calibration. Subject to the limitations of the speci- 


assure 


carbon 


stresses 


cobalt 


case, 


TABLE 2. S/ 


Specimen Pearlit 
Radiation } p/p 

, 
Co s! “4 + oo 1 14. 12> 
Fe ry 72 1. 30 7 
Cr 8 11 1. 20 14. 38 
Weighted averag« 12.74 
Theoretical 9.2 


men holder and Bragg angles, the angles of incidenes 
for the inclined exposures of each pair were 45°, 39° 
and 45°, respectively. 

The mechanics of getting the patterns, processing 
the films, reading and averaging the data, ete., wor. 
in no way significantly different with the steel speg. 
mens than with the The results, alone 
with the probable errors and the computed theoretieg) 
values, are presented in table 2. 

When examining this table, one should first notie, 
the variation in the observed probable errors. jy 
the stress; this error ts nearly three times as lare, 
for iron radiation for cobalt. The reason 
for this lies in the sizes of the Bragg angles givey 
in the table. The stresses measured by cobalt and 
by chromium radiation at the two larger angles 
however, do show the expected trend, which is the 
same as that in the theoretical values, in the 
order: pearlitic, spheroidized, and decarburized 
The mass absorption coefficients, u/p, of iron for 
the three radiations are also shown in the table. The 
cobalt radiation is the most penetrating, and js, 
therefore, the most reliable, because it reveals by 
diffraction the condition of the largest number of 
crystallites in the largest volume and is affected 
least by any residual surface preparation distortion 
that might be present. The average values of the 
stress shown in the table were weighted with respect 
both to the probable error in observation and the 
penetration of the radiation. 

The theoretical values in this table were computed 
with the aid of formulas for the specific lamellar 
and spheroidal cases given in La&szl6’s first paper 
The elastic constants used were those cited by Laszlé 
E—=30™ 10° Ib/in.2: v—0.30. the coefficient of linear 
thermal expansion of ferrite equals 14.98 107°/°C, 
and that of cementite 12.16 10°°/°C The tem- 
perature drop in the calculations was taken to be 
680° (. Again. the agreement between the theo- 
retical values and the observed averages Was reason- 


brasses. 


as it is 


ably good. 


4. Conclusions and Discussion 


In previous investigations involving X-ray stress 
measurement, it has generally assumed that 
carefully annealed material was strain-free. Laszl6é’s 
theoretical work showed that inequalities of thermal 
expansion would cause textural (tessellated) stresses 
in two-phase alloys, but no direct measurement 0! 


been 


@SSOS Lh steel 
Spheroidized Deearburiz Probabk 
error in 
stress 
7 o 
y} : yp lh/in 
TT 1.00« 10-4 3 56K 10 0.16«k10-4 0. 57K 10 1. 92X10 
1.59 6.84 2 45 10. 52 5. 34 
Is 10, 89 2.72 +, 29 2.02 
6. 1 58 2.53 
6.2 
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these stresses had been made. Contrary to the 
implication in one of LAszldé’s papers, the results of 
‘his work demonstrate that textural stresses of this 
type can be measured by changes in ‘the lattice 
constants if sufficient care is taken. The major 
experimental limitation is the problem of surface 
preparation, the routine of repeated cycles of light 
mechanical polishing and light electropolishing ap- 
parently being the best method available at present. 

In the present investigation, good agreement was 
obtained between theoretical stress values and those 
obtained by X-ray diffraction in the case of alpha- 
beta brass, and at least qualitative agreement in the 
ease of 1.05 percent carbon with various 
microstructures. The experimental precision was, 
in most cases, thought to be such that the probable 
i measurement was a few thousand 


steels 


error In stress 
pounds per square inch. 

In addition to the quantitative aspect of agreement 
with theory, these results reveal, by direct measure- 
ment, the impossibility of obtaining complete strain 
relief in a material that contains more than one 
phase. They will, perhaps, thus serve to emphasize 
an 
ever-present factor in stress corrosion, surface dis- 
tortion, fatigue, and similar problems. In_ this 
manner they may play some part in the very broad 
field of study relating the internal stress system in 
al metal specimen to the mechanical and physical 


the necessity of considering textural stresses as 


properties and behavior of the body. 


The authors wish to acknowledge the assistance 
of many with this work. 
{mong them are Hugh Logan, who supplied some 


persons hn connection 


of the early experimental brass allovs; personnel of 
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of 


the Naval Gun Factory and the foundry of 
Metallurgy Division at the National Bureau 
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timate brass material; Nathan Koenig, who did 
the greater part of the metallographie work; and 
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Dissociation of 4-Chloro-4’-Aminodiphenylsulfone 
Elizabeth E. Sager and Fleur C. Byers 


The ultraviolet absorption spectra of 4-chloro-4’-aminodiphenylsulfone has been de- 


termined. 


The dissociation may be considered as the acidic dissociation of the ammonium 


ion which is formed with the addition to the base of hydrogen ion from an acid such as hydro- 


chlorie acid. 


Spectrophotometric measurements of several concentrations of 4-chloro-4’- 


aminodiphenylsulfone at various known hydrogen-ion concentrations were made at 25° C, 


from which the values pA,, based upon the concentrations, were determined. 
pK, is 1.38, which corresponds to a value of 0.042 
The activity coefficient terms may be expressed simply as a linear 
Spectrophotometric measurements were also made at 15°, 20°, 


the value of pA, was derived. At 25° C, 
mole per liter for Ag. 


function of the ionic strength. 


25°. 30°, and 35° C on certain solutions in order to 
was found to be 21,000 joules deg-' mole 
24x 10-%, or pk , 1S 12.62 


1. Introduction 


Derivatives of diphenylsulfone, such as dihydroxy- 
diphenylsulfone and diaminodiphenylsulfone, have 
attracted attention as intermediates particularly in 
the synthesis of biocides. Heymann and _ Fieser 
have prepared many derivatives of the parent com- 
pound 4,4’-diaminodiphen) lsulfone in their reset rch 
studies of antimalarials [{1].'. They kindly supplied 
us with several compounds which they had prepared 
and we have used them in our studies of the spectral 
absorbance of organic compounds. The diphenyl- 
sulfones are almost insoluble in water, but solutions 
of the order of 107* to 10~°-.\7 may be obtained which 
are suitable for spectrophotometric measurements. 
Derivatives of diphenyvlsulfone absorb in the ultra- 
violet portion of the spectrum and measurements in 
this region afford a satisfactory means of observing 
some aspects of their chemical behavior. 

Compounds of similar molecular structure, specifi- 
eally benzophenone, 4-aminobenzophenone, and 4,4’- 
diaminobenzophenone have already been studied at 
the Bureau. Reports have been published including 
the ultraviolet spectra, dissociation reactions, and 
the effect of temperatures in the range of 10° to 40° C 
upon the spectra [2.3.4]. This paper deals specifi- 
eally with the acidic dissociation of the ammonium 
ion of 4-chloro-4’-aminodiphenylsulfone at 25° C, 
and with the effect of temperature upon its dissocia- 
tion, 


2. Experimental Procedure 


2.1. Materials 


t-Chloro-4’-aminodiphenylsulfone is a white crys- 
talline material that melts at 183 to 187°C [1]. The 
solubility in water is about 20 mg/liter at 25° C. 

Hydrochloric acid of reagent grade was used to 
adjust the hvdrogen-ion concentration of most of the 
solutions. Standardized sodium hydroxide, free 
from carbonate, was also emploved. Conductivity 
water was used to prepare the stock solutions and 
for subsequent dilutions 





From these, 


heat of dissociation. It 
constant, K,, is 


calculate the 


The basic dissociation 


2.2. Equipment 


A model DU Beckman spectrophotometer was 
equipped with a specially designed cell compartment 
with sides, bottom, and removable top which were 
water-jacketed. The temperature within the cell 
compartment was controlled by circulating water 
from a water bath controlled by a thermostat. 

A commercial glass-electrode assembly was used 
to make pH measurements of all solutions. 

All equipment except the water bath was housed 
in a constant-temperature room at 25°C. Solutions 
were prepared and stored in this room. 


3. Calculation of Dissociation Constants 
From Spectral Absorbance Data 


3.1. Absorbance Measurements 


t-Chloro-4’-aminodiphenysulfone exists as the base 
in water. When hydrogen ion from an acid such as 
hydrochloric acid is added, it is transformed to the 
ionic state, thus 
Cl SO. 


NH,+H* 


Cl SO, S—NH}. 
The reaction may be simply represented as the 
acidic dissociation of the ammonium ion, namely, 
RNHJ—RNH,+H"*. 

According to the law of absorption, at any given 
wavelength, 

Molar absorbance = A/bM (1) 
in which A is the absorbance of the sample being 
measured (—log,) transmittancy), 6 is the depth of 
solution in centimeters through which the radiant 
energy passes, and .V/ is the molar concentration of 
the absorbing compound. The quantity A/bM is 
frequently called molecular extinction coefficient, 
and is represented by the symbol «. Molar absorb- 
ance is expressed in liters mole~! em™', 
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FIGurReE l. Molar absorbance of 4-chloro- ,’-aminodiphe nylsul- 
fone in the molecular state RNH and the tonic state 


The values of molar absorbance for the molecular 
and the ionic states of 4-chloro-4’-aminodiphenyl- 
sulfone are shown in figure 1. The curve marked 
RNH, shows maximum molar absorbance of about 
17,000 liters mole! cm at a wavelength of 286 mu. 
The compound is transformed to the ionic state, 
RNH}?, in a solution of hydrochloric acid, 0.05-M], 
and exhibits maximum molar absorbance of about 
18,000 liters mole! em™' at a wavelength of 244 mu. 

Observations were made on four different series of 
solutions of 4-chloro-4’-aminodiphenylsulfone at 
various stages of the dissociation at known hydrogen- 
ion concentrations. A typical series is shown in 
figure 2. Absorbance values corresponding to 
experimental transmittancy values of a 3.75 * 107°-M/ 
sanidien of the diphenylsulfone in 1-cm absorption 
cells are plotted as a function of wavelength. Curve 
1 represents RNHz, in water, and curve 15, RNH; 
in 0.5-M hydrochloric acid. The intermediate 
stages in the dissociation are shown in curves 2 to 14 
inclusive. The two well-defined isosbestic pots at 
234.8 and 258 my not only indicate purity of the 
compound but show that only one reaction, that of 
dissociation, is taking place. The amounts of each 
species of the partially dissociated compound are 
proportional to their absorbance values (curves | 
and 15). The absorbance data may therefore be 
used to calculate the percentage dissociated into 
ions. Let a@ represent this percentage, which can 
be calculated from the following equation 


Agnu, A, 
a , ~ 
Apnn Apnn 


in which Apna, is the absorbance of the molecular 
species, A, is the absorbance at any dissociation 
step, and Agyn; is the absorbance of the ionic species. 








4 
_ 
' 
4 
‘\ —_ \\\ 
, XY 7 
i iX, 
Lt 1 —_ 
> A 260 280 aC 290 


WAVELENGTH, 17 

FIGURE 2. 
sulfone at 

iis dissociation (3.75 10 


Absorbance values of t-chloro-4 -aminodiphenyl- 
concentrations representing } 


5. V mn ] -« yn absorption ce Ils F 


various h ydroqge n-10n 


3.2. Calculation of Concentration Constants 


The concentration constant, A., can be calculated 
according to the equation 


kK [RNH,][H*} 
[RNH?] 
in Which [RNH,], [H*], and [RNH¢] are the concen- | 
trations of the free base, hydrogen ion, and the ionic 
species, respectively. Since a@ represents the degree 
of dissociation into ions, 1—a is the fraction of the 
remaining base, and we have the familiar equation 


1—a)/H | 


cx 


Kk 


It is advisable to observe the resulting a@ values 
for variations or trends in the calculated values. 
In this case it is obvious that best results can be 
obtained from about 280 to 292 mu. Calculated a 
values at nine w avelengths are given in table 1. 

Average values of @ are used to calculate each 
ratio of (l—a@)/a at the various hydrochloric acid 
concentrations as shown in table 2. The corrected 
hydrogen-ion concentration is that of the hydro- 
chloric acid less the amount of hydrogen ion used 
for the dissociation step corresponding to a. This 
figure multiplied by the ratio (l1—a)/a@ gives the 
concentration constant, A,. The negatives of the 
logarithms of the constants, or pA, values, are then 
found 


3.3. Estimation of pK, at 25° C 


The true thermodynamic constani, A,. is caleu- 


lated by eq (5 
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TABLE | {hsorbance values OJ chloro- '_aminodiphenylsulfone (3.747% 1075-M) and calculated values of a@ al several wave 
lengths 
Wavelength (my 
Soluti« A bsorbanct hl @ 
280) 282 2s4 285 286) 287 288 290 292 
fA 0. S884 0. 6056 0. 6162 0. 6108 0. 6216 0. 6216 0. 6198 0. 6091 0. 5035 
\a Mn 
{A 554 ASI7 5901 5U48 AW43 5043 AYLS 5850 A6R6 
2 la O44 043 046 045 04S O47 049 043 045 
fA 5339 A452 44 5575 975 575 5560 5AQS 5331 
la 104 110 109 109 112 111 111 105 110 
j {R28 {048 O45 5059 S072 \072 ROSY 4962 . 4828 
: a 201 201 197 199 144 10 149 200 201 
fA 4445 4541 109 4522 41 422 4609 4541 4413 
la 274 275 274 276 274 277 277 275 276 
fA $123 4202 4200 1283 4295 4295 4260 4202 1078 
la $35 337 336 335 335 334 338 $35 337 
i j830) 3893 3058 3969 3069 3069 3047 3803 3778 
la $41 393 38Y 390 392 391 392 390 392 
fA $575 3635 3075 3686 3686 3686 3675 3615 3507 
5 la 139 140 130 130 441 440) 440) 439 441 
{A $372 3410 3458 3468 3468 3468 3449 $391 3208 
a 17s 4s1 i7s 178 179 i7s 479 17s 479 
i 3179 $224 {261 3270 3270 3270 3251 3188 3098 
la 5} 514 12 512 513 512 514 514 515 
(4 OR48 yt] 2007 2016 2916 2007 RO) 2840 2757 
la 7s SSO 575 574 575 57 a7 576 577 
i 2388 23H 2411 2411 2411 2403 23585 2358 2284 
- a it tt Hos Ho3 663 663 Ho4 H6H2 663 
| 1931 1904 1904 1891 1891 Iss 1878 Is44 1791 
la 752 754 752 753 754 74 753 753 753 
. fl 1209 1163 113! 1124 1118 1107 1096 1085 1046 
. a Su) SAU SAS SAS SAS SSU KU SAT RAS 
i O630 O50 O50 0482 0477 0467 (462 0448 0429 
a 1. 000 
TABLE 2 Dissociation of 4 chloro- ’_aminodiphen ylsulfone (3.747 > 10 5. M) 
Corrected 
: AV l xia Molar concen concentration Concentration logioA loz f=1.0y4 | pK. pk. 
I RVU tration of HC] of hydrogen constant, A pk log f) 
ion |H* 
wx s 
2 ne m +) 2 wz i 2 OOO 10 4 148x110 1. 382 0, OO2 1. 380 1. 005 
; ws 174 5 Os 5 OO1 4+ OSS 1, 3&8 O05 1. 383 1.012 
j aM + (Wn) 10. 010 10. 008 + OO 1. 308 O10 1, 388 1.024 
27t » 2 15.015 15. O05 , O36 1. 405 O15 1. 380 1.035 
2). O20 A). OOF 3. 953 1. 403 020 1. 383 1.048 
) F 25. O25 25. 010 4, 804 1.410 . 025 1. 385 1. 059 
8 44 1.2 30. O30 40. 014 3, 820 1. 418 . 030 1, 38S 1. 072 
) 174 (is 35. O35 35 O17 , 806 1.420 035 1. 385 1.084 
449 40. 040 1), 021 3. 798 1. 421 040 1. 381 1. 097 
50. O50 50. 028 3. 682 1. 434 050 1. 384 1. 109 
2 Ons s 70. 000 69. O75 3. BAS 1. 449 070 1.379 1.175 
sos 100. 000 a9. O72 3. 279 1. 484 100 1. 384 1. 259 
$ xSS 200). OOO 109, O67 2. 540 1. 505 2) 1. 395 1. 585 
wn c SOO. OOO 
. Justeen When the values of pA, are plotted as a function of 
A Kk, a the ionic strength, w, a satisfactory extrapolation to 
RNH 


in Which the activity coefhicients are identified with 


appropriate subscripts For simplicity this com- 
bined term will henceforth be referred to as f. Then 
pk pk log f 0 


The slope of the 
Therefore 


the value of pK, can be made. 
line is about 1.0, as shown in figure 3. 

log f may be considered equal to 1.04. The values 
of log f may be subtracted from the experimental 
pK, values to obtain individual pK, values that can 
be averaged. By this method of calculation pK, is 
found to be 1.38, with an estimated uncertainty of 
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Figure 3 pk, values plotted as a function of the ionic strenqth Ficure 4 ph values at temperatures from 15° to { 
for extrapolation lo pk, plotted as a function of 1/T 
TABLE 3 Calculation of pk, from absorbance data at 15° to 35° ¢ 7 
pro} 
he 
rem 
pera Absorbance at 286 ing O.015-M Ht 0.035-M HC] ) and 
ture allio 
only 
( 1RNA innuy 1 onsen 1 osHci a)a A pk pRK« l—a)ia A A pK of | 
15 ©. 5800 0. 0366 0. 3968 0. 2831 7 204 1. 532 1.51 0. S2S » 2 1. 53s 1. we 
a” 5703 O410 4144 {01 24 033 1.472 1. 457 2s 032 {XS 1. 453 tem 
25 5767 0453 $300 $224 2. 623 303 1. 406 1. 391 1.000) 0382 417 1. 382 T 
3 725 0520 14 $400 O44 M42 l m0) l 4 uy (43 1, 328 } 
$5 w4 O575 $522 ar is7 od 281 att Ht 1 "T Is U 
inst 
liqu 
— , _ , ' - — the 
+0.01 at 25° C. Values of f corresponding to each | the pA, values against 1/7>< 10°, as shown in figure 4, bra 
value of the ionic strength in this series are shown in | vields a practically straight line indicating that the in | 
the last column of the table heat of dissociation is nearly constant over the tem- the 
perature range studied. A/7/ is found to be approxi- ma 
3.4. Heat of Dissociation mately 21,000 joules deg™! mole~'. This value is not tere 
; ah very different from one of 19,000 reported earlier in leas 
he heat of dissociation is calculated from the | a study of 4-aminobenzophenone which was based me’ 
equation on more extensive data. pla 
d log K/dT=All/2.303RT- 7 val 
: , ;, 3.5. Basic Ionization Constant, K ran 
in which A// represents the change of molar heat =. cal 
content for the dissociation, A is the gas constant ry - “t . ys . , hs 
ys : , 1 -1 a The basic ionization constant, A>, is derived from a 
(8.3144 joules deg! mole~'), and 7 is the tempera- .s ee 
. . a “Se ee 979 12 the usual relationship ani 
ture in degrees Kelvin (degrees Celsius plus 273.16). rep 
Four solutions in one series, 3.48 10™°-M/ with , > t] 
Pa. , K,= Kh, Kg, 6 . 
respect to 4-chloro-4’-aminodiphenylsulfone, were 
observed throughout the ultraviolet spectral range — . . r _ , ’ 
— ; 20 4 on° 9 in which A, is the activity product of water. For 4- 
at five temperatures, namely 15°, 20°, 25°, 30°, and , . . “1 ys he 
ae a | ae chloro-4’-aminodiphenylsulfone Ay is 107'/4.2 10" 
35° ¢ They represented the dissociated state, the or ’ 9 “1 pn . : . 
* Ae Ar te, . . - | at 25° C, or 2.410-"%. This gives an equivalent . 
undissociated state, and two intermediate stages of value of 12.62 for vk ral 
- pipe: hn = . alue 2.62 I IIA». 
dissociation at 0.015-M and 0.035-M hydrochloric / th 
acid concentrations. Absorbance values at the dif- du 
ferent temperatures, at wavelength 286 my only, are 4. References eT) 
shown in table 3. The K, and pK, values are cal- | bk 
culated for the two concentrations in the usual man- | [1] Hans Heymann and Louis F. Fieser, J. Am. Chem. So th 
ner. In this limited temperature range no significant Sande 1979 (1945 rm) re 
Tae iy RPT : “ eer : [2] Elizabeth E. Sager and Iris J. Siewers, J. Research NBS 
error is introduced by considering concentrations in- 45. 489 (1950) RP2162 m 
dependent of temperature. Assuming that the same | [3] Elizabeth E. Sager and Iris J. Siewers, J. Research NBS pu 
activity coefficient terms (—log f=1.0 wp) apply at 49, 7 (1952) RP2337 IS 
the temperature used, which is, of course, not strictly | 4) ey r at A Fleur C. Byers, J. Research NBS fre 
“ , : . : a8, 33 O57 . P27: t 
true, pA, values can be estimated from data for the oe > lat 
two solutions. They agree within 0.01 pA unit at th 
each temperature, so an average is used. A plot of WaAsHINGTON, March 19, 1957. ps 
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Speed of Sound in Water by a Direct Method’ 


Martin Greenspan and Carroll E. Tschiegg 


The speed of sound in distilled water was measured over the temperature range 0° to 


100° © 
nomial and in tables 


with an accuraey of 1 part in 30,000. 
The water was contained in a cylindrical tank of fixed length, termi- 


The results are given as a fifth-degree poly- 


nated at each end by a plane transducer, and the end-to-end time of flight of a pulse of sound 
was determined from a measurement of the pulse-repetition frequency required to set the 


successive echoes into time coincidence 


1. Introduction 


The speed of sound in water, ce, is a physical 
property of fundamental interest ; it, together with 
the density _ determines the adiabatic compressibility, 
and eventually the ratio of specific heats. The vari- 
ation With temperature is anomalous; water is the 
only pure liquid for which it is known that the speed 
of sound does not monotonically with 
temperature. 

There is also a practical interest in ¢ in that water 
is used as a standard liquid for the calibration of 
instruments that measure the speed of sound in 
liquids automatically, both in the laboratory and in 
the field. In fact, it was in connection with the eali- 
bration of such ‘‘velocimeters” [1]? that our interest 
in this work was first aroused. In the first place, 
the available data scatter widely, as recent sum- 
maries [2, 3] clearly show. In many cases, the dis- 
erepancies far exceed the claimed accuracy or at 
least the precision of the methods, even when the 
methods compared are the same. In the second 
place, there exists no set of data that gives a smooth 
variation with temperature over any considerable 
range. In particular, the best of these data vield 
calibration curves for our velocimeters which are 
badly curved instead of straight (as they should be), 
and about which the data scatter irregularly, but 
reproducibly The results here presented are free of 
these object ions. 


decrease 


2. Method 


At the top of figure 1 is a schematic of the appa- 
ratus. The sample is confined in a tube of which 
the ends are plane, parallel, electroacoustic trans- 
ducers, quartz crystals in this case. If the left-hand 
crystal, say, is excited by a short pulse from the 
blockmg oscillator, the oscilloscope, which measures 
the voltage on tne right-hand crystal, will show a 
received pulse and a series of echoes, as indicated 
in idealized form on the line below (fig. 1). The 
pulse repetition frequency of the blocking oscillator 
is controlled by a sine-wave oscillator, and if this 
frequency were adjusted so that each blocking oscil- 
lator pulse coincided with the first received pulse of 
the next preceding evele, then the oscillator period 
would equal the time of flight of the pulse. However, 


This work was supported in part t the Office of Naval Research under 


contract N A-onr-70-48 


Figures in brackets indicate the literature references at the end of this paper 


as the two pulses have different shapes, the accuracy 
with which the coincidence could be set would be 
very poor. Instead, the oscillator is run at about 
half this frequency and the coincidence to be set is 
that among the first received pulses corresponding to 
a particular electrical pulse, the first echo correspond- 
ing to the electrical pulse next preceding, and so on. 
Figure 1 illustrates the successive signals correspond- 
ing to three electrical input pulses. The input pulses 
fall halfway between the a for which the coinci- 
dence is set, so that they do not tend to overload 
the amplifier or distort the oscilloscope traces. The 
period of the oscillator, when properly set, multi- 
plied by twice the length of the tank, is the speed of 
sound in the sample. 

The oscilloscope trace actually looks like that 
shown in the inset (fig. 1). The first cyele corre- 
sponds to sound reflected from the inner faces only 
of the transducers, whereas the succeeding cycles 
correspond to sound reflected one or more times from 
an outer face. Therefore, the coincidence is set by 
maximizing the peak on either the first or second 
balf-eycle; the same result is obtained in either case 
but the second half cycle is easier to use because it 
is bigger. What we are measuring here is the speed 
corresponding to the first arrival of the signal; in a 
nondispersive liquid this is the same as the phase 
velocity. It is true that the coincidence is made at a 
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FIGURE 1. 





Schematic of method. 


The three lower lines show in idealized form the events corresponding to three 
successive electrical pulses Che short, thick line represents the input pulse 
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time one-fourth or three-quarters of the transducer 
period later than the time of first arrival, by which 
time there is opportunity for sound traveling by 
paths other than the shortest to affect the location 
of the maximum. However, the results are inde- 
pendent of whether the first or second half cycle is 
used; they are also not affected by substituting 
cry stals of twice the thickness, or by changing the 
diameters of the tank, or of the hot electrodes. These 
results lead us to believe that the error introduced by 
this maximization technique is negligible. 

The question has been examined also in another 
way. Suppose a coincidence to have been made at 
frequency f; others can then be made at submultiples 
of f. At the frequency //2 for instance, the first 
received pulse corresponding to a particular input 
pulse coincides with the second echo (not the first, 
as before) corresponding to the electrical pulse next 
preceding, and so on. Effectively, the sound pulse 
is timed over a path twice as long as before. It is 
found that the measurements at f and near //2 are 
substantially identical, so that the error in question 
is less than, or at most comparable to, the experi- 
mental error of the time measurement. 


3. Apparatus 
3.1. The Delay Line 


The disassembled delay line is shown in the photo- 
graph, figure 2. The length of the tank is about 
200 mm, and the bore about 13 mm. The filling 
holes are sealed by plugs having Teflon gaskets; a 
small hole in one plug provides pressure release. 

The tank is of a chromium steel ® which, after heat 
treatment, takes a good optical finish. Because 
this steel is not so corrosion resistant as the nickel- 
chromium stainless steels, the bore of the tank was 
heavily gold plated. 

The ends of the tank are optically flat and parallel 
to within less than 1 yw. To these ends are carefully 
wrung the 0.8-mm thick x-cut quartz crystals, which 
also are optically flat. The caps, when bolted on, 
clamp the crystals through neoprene O-rings. A 
coaxial cable passes through a seal in each cap, and 
the center conductor makes contact with the outer 
(hot) electrode of the crystal through a light spring. 

The outer electrode is a 9 mm circle of aluminum- 
backed pressure-sensitive adhesive tape. The inner 
(ground) electrode is of fired-on gold and is about 
12 mm in diameter. Contact is made through a 
light gold-plated helical spring which touches the 
electrode around the edge and bears on & shoulder 
machined into the bore. The inner electrodes and 
springs are unnecessary if the sample has high con- 
ductivity or a high dielectric constant; they are 
usually omitted for water and aqueous solutions of 
salts. Figure 3 is a schematic drawing of one end 
of the assembly. 

_ The length of the tank was measured at 20° C 
and the coefficient of thermal expansion of the steel 
was measured on a sample cut from the same bar as 


eel 


Firth-Sterling type}B-440A 





Delay line, 


Above the tank are the plugs which close the filling holes. 
the caps through which pass 
crystals seen in the foreground 


Fiat RE 2 disasse mbled. 
snd at the ends are 


the electrical cables and which alse Clamp t 



































Schematic of one end of the lank, 


FIGURE 3 showing the crysta 


and cap Se mbly. 
1, Tank; 2, cap; 3, plug sealed with 


through pressure relief tulx 1, quartz cr 
make contact to electrode 5 neoprene 


reflon O-ring (thermocouple 
ystal wrung onto end of tank (springs 
O-ring; and 6, coaxial cable. 


Passes 


was the tank and heat-treated together with it 
From these data, the length of the sound path Is 
known to better than 2 parts in 10° at any tempera- 
ture between O° and 100° C. It is, of course, 
necessary that the crystals be w rung down with 
great care so that the fringes disappear all around 
the periphery, to achieve this accuracy. The clamp- 
ing gaskets must bear directly over the contacting 
surface and not spread out over the unsupported 
area, else the crystal will bend. With these pre- 
cautions, the delay line may be disassembled and 
reassembled repeatedly with reproducible results. 
If the crystals have been properly wrung on and 
clamped, they cannot be removed by hand after 
several days, but must be soaked off 
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3.2. The Electronics 


The electronic circuits (fig. 1) are, for the most 
art, conventional. However, the oscillator and the 
forming circuits must be exceptionally free of 
jitter. In addition, the oscillator must be provided 
vith fine frequency control, so that it can be set 
within the required sensitivity of measurement, and 
it must be so stable that the lrequency does not 
change during the counting time by enough to alter 
the count by more than one. 

The blocking oscillator produces a pulse about L100 
y high and 0.05 to 0.25 ysee wide. It is best driven 
by a large, fast pulse such as is gotten by differentia- 
tion of a square wave derived, in turn, from the sine- 
wave generator.* The jitter may be reduced further 
by means of a narrow band filter after the oscillator. 
“The receiving circuit consists simply of a short 
length of low-capacitance cable, a wide band (5.5 
\[e. in this case) amplifier of gain 100 to 1,000, and 
» high-frequency type oscilloscope equipped with 
fast sweeps. The sweep is triggered from the 
oscillator through a variable delay; the necessary 
delay time is about half the oscillator period. 


pulse- 


3.3 Temperature Control and Measurement 


The delay line, suspended from its cables, is 
deeply immersed in a 27 gal, well-insulated, water 
bath. The bath is provided with 2 pump-type 
stirrers, 3 heating coils, and a cooling coil connected 
to a small refrigeration unit. The temperature of 
the water adjusts itself so that the losses equal the 
power input to the heating coils, and various temper- 
atures are obtained simply by varying the power input. 
The temperature is, by this means, easily held to 
within less than 0.005 deg C for the interval of time 
required for the measurements, except that above 
75° C, or the variation may 
0.03° C; at the higher temperatures, however, the 
thermal coefficient of the speed of sound in water is 
rather low. 

The temperature of the bath water is measured 
with a platinum resistance thermometer and Mueller 
Bridge. A differential thermocouple has one junc- 
tion in the sample in the delay line, and one tied to 
the platinum thermometer; it passes through the 
pressure release tube (fig. 2). The thermocouple 
reading serves to indicate when the sample and bath 
Water are substantially in thermal equilibrium, and 
measurements are made when the discrepancy is less 
than 0.01° C (somewhat greater at the high temper- 
atures). The thermocouple and galvanometer com- 
bination was calibrated for small temperature dif- 
ferences against the platinum resistance thermom- 
eter so that small corrections to the temperature 
readings could be made. 


sO, 


3.4. Samples 


The measurements here reported were made on 
three separate samples of water. One sample was 
ordinary laboratory distilled water. This was boiled 


ED 
‘ The square-wave generator must be of the type 


that amplifies and then clips 
the input sine wave he free-running ( 


i type is not suitabk 


become 0.02° or 


and poured, while still hot, into the preheated tank. 
Although dissolved air has a negligible effect on the 
speed of sound in water [4], it is desirable to exclude 
air and so prevent possible bubble formation on the 
transducers. 

The other two samples were vacuum distilled di- 
rectly into the tank. The tank was placed in an ice 
bath and connected to a flask of distilled water. 
The system was then evacuated, and the water 
allowed to distill over at about 50° C. 

The results of the three runs were the same within 
the errors of measurement; the data were, therefore, 
combined 


3.5. Technique 


The water bath was cooled to just above 0° C, 
and the heaters were operated at low power to sta- 
bilize the temperature. (Below room temperature 
the refrigeration machine was run continuously.) 
After the readings were taken, the power input to 
the heaters was increased, and so on until the tem- 
perature was just below 100° C. When the tem- 
perature was stabilized, as indicated by the con- 
stancy of the Mueller bridge reading and the near 
zero reading of the thermocouple galvanometer, the 
coincidence was set on the oscilloscope by one ob- 
server and the frequency (doubled for convenience) 
was measured by counting cycles for 10 sec (about 
75,000 counts) by means of an electronic counter. 
At the same time, another observer balanced and 
read the Mueller bridge and read the thermocouple 
galvanometer deflection. 

While the temperature readings were being made, 
the coincidence was independently set and the re- 
sulting frequency measured three times or more. 
The various readings were always within the +1 
count inherent error (even for different observers) 
and the modal value was recorded. This, divided 
by 20 and multiplied by the length of the tank at 
the particular temperature, was taken as the speed 
of sound, ¢, corresponding to the temperature, 7’, 
obtained by calculation from the platinum: ther- 
mometer and the thermocouple readings and the as- 
sociated calibration data. All temperature calcula- 
tions were made to the nearest 0.001° C and the final 
result was rounded off to the nearest 0.01° C. 

In order to insure that the coincidence was set on 
the proper cycle, it was first set approximately, using 
the coarse frequency control, at a moderate sweep 
speed and low oscilloscope gain, so that the entire 
pulse was visible on the screen. The sweep speed 
was then increased while the delay was readjusted 
to keep the proper cycle centered. Next, the gain 
was increased while the base line was moved off the 
screen to keep the point of extreme deflection cen- 
tered, and the amplitude was then adjusted to a 
maximum using the fine frequency control. 


4. Results 


From readings taken at 83 temperatures between 
0.14° and 99.06° C, the calculated values of the 
speed of sound were fitted by the electronic com- 
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FIGURE 4 Deviations, r, of equation 1 from the data 
puter SEAC by the method least squares, to a significant figures, so that on account of rounding-of 
fifth-degree polynomial, errors, the tabulated differences in some cases differ 
5 by one unit in the last decimal place from the dif. 
c=) 0,1". | ferences of the tabulated values of ¢. It is believed 
=e (see section 5) that the systematic errors do not 
. ° . . ’ ay ay Tr j "2 7 ’ ‘ Ss < _ 
lhe reduction in the residual sum of squares over a a ? oe) ee ine folk Phe tabl — there. 
. . . _ vi » se » r ‘ i ¢ . 
fourth-degree polynomial, due to fitting the fifth- | 10": D¢ Used In the folowing manner. In table 2 
: linear interpolation should be performed to th 


degree term, was statistically significant at a prob- 
ability level less than 0.005, and the deviations of 
the data from the fifth-degree polynomial showed no 
statistically significant indication of lack of random- 
The deviations are plotted against tempera- 
ture in figure 4. 

The values of a; in eq (1 


hess. 


), for ¢ in meters per second 


(m/s), and 7 in degrees C, are: a,=1,402.736; a 
5.03358; ds 0.0579506; ag=3.31636 & 1074: a, 
1.45262« 107°: and a,=—3.0449*10°°. The stand- 


ard deviation of the measurements is 0.0263 m/s, or 


about 17 ppm. Estimated standard deviations of 
the values of ¢ predicted by eq (1) were calculated 
for five representative temperatures. The results 
are given in table 
TABLE 1. Estimated standard deviation (s. d.) of values o 
predicted by equation ] 
Pemperature 
( mis ppm 
0.0114 81 
10 OOS 4. § 
5 “ ) s 
70 ONL 10 
100 O45 4 
Table 1 and figure 4 make it clear that eq (1 


Ssatis- 
intro- 


together with the listed constants, provides a 
factors interpolation formula, and the errors 


duced by its use are small relative to the possible 
svstematic errors of measurement (see section 5 
The values given in tables 2 and 3 were calculated 


from eq (1) by SEAC 

Table 2 vives the speed sound mn meters pel 
second for each degree C from 0 to 100, and table 3 
gives the speed of sound in feet per second at inter- 
vals of 2 deg F from 32° 212°. In each case, the 
differences, which are listed for convenience in inter- 
polation, were calculated from a table having more 


to 


nearest 0.01 m/s and the final result rounded off to 
the nearest 0.1 m/s. The error will then not exceed 
one-half unit in the last place, i. e., 0.05 m/s. Linear 
interpolation in table 3 will vield errors that do not 
exceed 2 units (0.2 fps) in the last place. 


5. Discussion 


Following is a list of the known possible sources 
of error and an estimate of the upper limit of each 
error 

5.1. Frequency 


As already stated, the frequency was measured by 
counting cycles for 10 sec; the total count was 
about 75,000. The inherent error is +1 count, but 
in all cases the mode of at least three independent 
readings, of which, at worst, two were the same and 


the third different by one, Was taken as the observed 
value. The counting error can thus be as great as 


1 part in 75,000, but as it is random, the effect on 
the final results is negligible, as indicated in section 4 
The 10-see time base was obtained by division from 
a 1-Me crystal oscillator which is stable to 2 parts 1D 
10’ per week, and which was compared with signals 
from WWYV or from a local precision sti andard. The 
errors due to inaccuracies in the time are, 
therefore, also negligible. 


5.2. Length of Path 


The length of the tank across its polished ends at 
20° C was determined within I, i. e., 5 ppm. Ther- 
mal expansion measurements were made at 20°, 60°, 
and 100° C; the lengths at intermediate tempera- 
tures were calculated by quadratic interpolation 
The maximum absolute error in the thermal expan- 
sion coefficient is estimated at 0.2 ppm; this acecu- 
mulates to 4 ppm at 0° C, and to 16 ppm at 100°C 


base 
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i TABLE 2. Speed of sound in water, metric unis 
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Thus. the total uncertainty in the length of the | attention to avoidance of clamping pressure too near 
tank is about 5 ppm at 20° C, and increases with | the unsupported areas of the erystals, the crystals 
1. al temperature both Ways; al 0° ( it becomes about might deflect enough to cause very sizable errors. 
oe 9 ppm and at 100° C, about 21 ppm. The present design makes it possible to disassemble 
Ci- Tt . . ° 
e0° Phe question arises as to how closely the length | and reassemble the delay line repeatedly without 
of, ol the sound path inh the sample, . Oi. the distance affecting the result by a detectable amount; this 
era- | : ‘ : . : 
ow between the inner faces of the transducers, approxt- holds true when the crystals normally used, which 
pan- mates to the length of the tank across the ends to are 0.8 mm thick. are replaced by crystals 1.6 mm 
ae which the ervstals are wrung. Experience with | thick. It, therefore, appears that errors produced 
i- P ° . ° 
o(" developmental models showed that unless the assem- | by misplacement or deformation of the crystals are 
bly were very carefully made, with particular insignificant 
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5.3. Setting the Coincidence 


As explained in section 2, it is believed that no 
measurable error is introduced by the technique of 
maximization of the second half evele of the received 
pulse. However, a word should be said about the 
effect of personal bias on the part of the operator. 
The operators report, to varying degrees, tendencies 
to adjust not only for maximum height of peak, but 
also for maximum symmetry and sharpness of peak. 
Long experiment has convinced us that any of the 
three criteria lead to sensibly the same result, so 
that although different operators weigh the three 
criteria differently, they reproduce each other’s 
settings so well that the discrepancies are negligible 
relative to other sources of error. The assumption 
is implicit that the errors of bias do not much exceed 
the discrepancies among individuals. 


5.4. Temperature 


The Mueller bridge with which the resistance of 
the platinum thermometer was measured has a least 
count of 0.0001 ohm corresponding, for a 25-ohm 
thermometer, to about 0.001° C. The bridge was 
calibrated internally so that the indicated resistance 
in terms of the internal standard is correct to about 
0.0002 ohm aside from temperature effects and slow 
drifts in the arm ratio and in the zero. Allowing 
for these, it is estimated that the bridge error does 
not exceed 0.005° C; errors in the calibration of the 
platinum thermometer itself and those due to heating 
by the bridge current are much smaller. More im- 
portant is the error that arises from thermal gradients 
in the bath. On the assumption that this does not 
exceed half the reading of the differential thermo- 
couple which, it will be recalled, measures the dif- 
ference between the temperature of the platinum 
element and that of the sample, an upper limit to the 
corresponding uncertainty in the speed of sound, 
ec, was calculated at various temperatures from the 
known thermal coefficient of ¢. This upper limit 
is zero at 74° C, where ¢ is stationary and increases 
steadily in both directions, reaching about 25 ppm 
at 0° C, and about 14 ppm at 100° C 


5.5. Purity of Sample 


Because the results obtained on ordinary labora- 
tory distilled water were indistinguishable from those 
obtained on the same water redistilled in vacuum 
directly into the apparatus, it Is felt that the remain- 
ing impurities do not have a measurable effect 
Several measurements made on local tap water gave 
results about 30 ppm higher than for distilled water 


5.6. Over-all Accuracy 


From the foregoing discussion it appears that the 
major sources of error are the uncertainties in the 
length of the path and in the temperature. Both of 
these are temperature dependent; their sum is an 
upper limit to the total error. This is about 35 
ppm at 0° C; it falls to 15 ppm at 40° C and is almost 


constant at this value out to 70° C, and rises ty 
about 35 ppm at 100° C It is upon these considera. 
tions that the recommendations for the use of thy 
tables in section 4 are based. 

The values of ¢ here reported are lower than thos 
of most other workers, in particular the value a 
30° C is about 0.4 m/s below that of Del Grosso 
Smura, and Fougere [3], whose work with the ultra. 
sonic interferometer is perhaps the most carefully 
planned, executed, and analyzed work of this typ 
to date. It was, therefore, felt desirable to perform 
an independent experiment using an apparatus and 
a method as different as possible from those of bot} 
Del al., and ourselves. An apparatus 
was constructed with which it is possible to measure 
as a function of distance, the phase on the axis of g 
beam of progressive waves emitted by a small piston- 
like radiator. If the wave were plane, the phase ¢ 
would vary linearly with distance «, and the phase 
speed c would be 2afa ta) where f is frequeney, In 
the present Case, the wave Is not plane and the slope 
of the curve 22fr versus ¢ depends on x and on the 
reometry of the arrangement. However, the theory 
enables us to select a distance sz) of the receiver from 
the that the departure of 
2afdrde from ¢ is as small as desired 

Five runs were made in distilled water at tempera- 
tures between 15 4 The prircipal un- 


Grosso, et 


source such for r>v, 


and 25° © 
certainties are thought to be first, one of about 40 
ppm corresponding toa possible error of O.OL° C in 
the temperature, and second, one of about 56 ppm 
related directly to the inaccuracies of the screw wit! 
which the receiver displacement was measured 
These are independent. However, in the worst Cast 
of the 5, the result differed from the value gotten 
from table 2 by only 27 ppm. The value of Del 
Grosso, et al. [3] disagrees with that of table 2 by 272 
ppm 
This work will be reported in detail elsewhere 


The authors are grateful to the personnel of the 
Engineering Metrology Section and of the Length 
Section, in whose Laboratories the length of the tank 
and its thermal expansion, respectively, were meas- 
ured. Thanks are particularly due to Joseph M 
Cameron of the Statistical Engineering Section who 
advised the authors on problems of data processing, 
and who performed the curve-fitting computations 
on SEAC. 
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Pattern of a Flush-Mounted Microwave Antenna 
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Wait 


The numerical results for the far zone radiation from an axial slot on a circular cylinder 


of perfect conductivity and infinite length are 
large diameter evlinders can be expressed in a 
enaleulations for arrays of slots on a gently 
related diffraction problem considered by Fock, 


1. Introduction 
The radiation from flush-mounted microwave 
antennas has aroused considerable interest in recent 
vears [1—6}.’ Krom a_ theoretical standpoint the 
problem would seem to be simple enough; the surface 
on which the antenna is mounted is considered to be 
adequately represented by a sphere or cylinder 
whose radius of curvature is matched to that of the 
local curvature of the actual surface. Unfortunately, 
the formal solution of the problem, in terms of the 
classical harmonic series involving integral or half- 
integral order Bessel functions, always converges 
very poorly when the radius of curvature is large 
compared to the wavelength. There are two well- 
known alternative representations |6, 7|, however, 
which can be used to advantage In certain 
In the illuminated region of the antenna, geometrical 
optics is most satisfactory, while deep in the shadow 
the rigorous harmonic series can be converted into 


cases. 


the rapidly converging residue series. 

The calculation of the radiation field of a flush- 
mounted antenna in the tangent plane (the classical 
light-shadow boundary) is not readily treated by 
either geometrical optics or the residue series. In the 
former case the field is indeterminant, and in the 
latter case the convergence is extremely poor and 
would actually diverge in the illuminated region. 
Despite the fact that the harmonic series is cumber- 
some, it is valid in this transition zone between the 
illuminated and shadow There- 
fore, it is desirable to attempt to adapt the harmonic- 
series representation to surfaces of large radius of 
curvature. This is the purpose of the present paper. 


regions of space. 


2. Theoretical Basis 


A thin axial slot, cut on a circular cylinder of 
infinite length and perfect conductivity is con- 
sidered because it gives rise to a plane-polarized 
radiation field. The eylinder is taken to be of radius 
a and coaxial with a cylindrical coordinate system 
00,2). The slot that extends from 2, to 2. at ¢=0, 
has a voltage distribution V throughout its 
length. The radiation pattern is best expressed in 
terms of spherical coordinates (7.6.6) where @=0 


~) 
{ 


t the end of this paper 


255 


curved surface 


discussed. It is shown that the results for 
universal form that is suitable for pattern 
The work is compared with a 


corresponds to the axis of the cylinder. It has been 
shown [9] that the electric field, which has only a @ 
component, is given by 


iar 
’ é ’ 
EF, S(0)M(x.¢), (1) 
p 
where 
k sin 6 (”., ; 
S(0)= Viz)e** 2? dz, (2) 
2r 2 
1 =. «.e°'"'* cos me 
M (x,¢) >> TL ? (3) 
WL mao H'?)"(z) 


r=ka sin 0, k=2r/free space wavelength, e)=1,€, 
2(m 0), and 72’ (x) is the derivative of the Hankel 
function appropriate for a time factor exp (io). 
Equation (1) is valid for kr sin @>1. 

The function S(@) is the space factor of the slot, 
or an array of collinear slots if they were cut in an 
infinite plane conducting sheet. The function 
M(z,¢) is called the “cylinder space factor’ as it 
fully describes the effect of the finite diameter of 
the eylinder on the radiation pattern in both the 
@ and @ directions. For purposes of computation, 
it is written 

Mi(x.0) M(x,0) fa? (4) 

The phase function a(z,¢) is a rapidly varying fune- 
tion of @ for the larger x values. For this reason, it 
is desirable to express a in the illuminated region 
(0< @ <r/2) as a geometrical-optical term plus a 
correction factor A(z,@), as follows [9] 

a(r,o@)=a(zxr,0) —xr(1 —cos ¢) + A(z,¢). (5a) 
In the shadow region (2/2<\¢\<), it is desirable to 
express a as a physical-optical term plus a correction, 
which is also designated A(z,¢), in the following 
manner. 

a(zr,o)=a(r,mr/2) +-(27/2—)2r+Alz,¢). (5b) 
The amplitude |\/(2,¢)| and the phase correction 
A(x.) of the cylinder space factor are shown plotted 
in figures 1 and 2 for @ from 0 to 180° for x=10, 12, 
15,18,and21. To prevent troublesome overlapping, 
the ordinates are shifted for each curve by a constant 
amount. 





i ST 2 Gs GRE SG ee GA GR For @<90° but not near 90°, the value of the 
™ space factor is quite well approximated by the 
\. ! geometrical-optics approximation. That is 


< A(r,.o) ~0 
and 
S so ’ ] Mir.) ~ | 


1+— oe Deep in the shadow, that is ¢) >90° but not near 
90°, \M(2,¢) has the form of a damped standing waye. 
Te and can be characterized by a function of the form 


u \ \ 7 


> \~U where 
“ nh A(p ~exp | it(wa rd) B 
\ with 7=0.808 exp(—/m/3). This exponential form 
for .1(@) is the first term of a rather complicated 
] residue series {7 It is only valid for ka>1. The 
: -. higher order terms (higher modes) are also parametri 
\ ey in (wa/dX)*8, where B is some angular distance. It 
~. is, therefore, logical to expect M I.d) to be parametric 
YA / in (2/2)'8(@—-w/2) and (2/2)[(32/2)—o)|. In facet. if 
the ripples are imnored., the function \/ I .@ should 
be a function only of (2/2 o nr) 2 The dotted 
_ curves in figures l and 2 are believed to be the 
appropriate form of \/(7.¢) when the standing-wave 
effect, due to the interfering traveling wave from the 
other side of the evlinder, is removed. These Fy 
“smoothed” curves are then replotted as a function 
of .\, where JA\ r/2 o—n/2 The results are 
shown in figure 3a where both the values of amplitude 
MCX and the phase correction ACA) for 2=10, 
iz, ia BO and 21 fall on the same set of curves for 
the range of .\V indicated The amplitude is also 
shown In figure 3b. being plotted ona log scale. 
Having the space factor for the evlinder plotted 
4 -™ in this universal form suggests that one may ex- 
. \ trapolate the results to larger values of «. This 
r concept, embodied in presenting the results in terms 
= j j of a parameter XY, is not unrelated to the notion of 
a . ey , angular distance that has been successfully emploved 
+~f\° in the representation of field strengths of a trans- 
- ~ | mitter at large ranges over a spherical earth 10 
o a. tT In the latter instance, the angular distance 1s de- 
ai \ iy fined as the angle subtended at the center of the 
- | earth between the horizons of transmitting and 
° “ U receiving antennas. In the present situation, the 
r mn f} angular distance is @—7/2 since the source antenna 
AY} | is on the evlinder and the observer is at infinity. 
™~ J ry While the present problem is concerned directly 
: with the radiation of an axial slot on a cireular eylin- 
VA} /P der, there is a direct application of the results to the 
reciprocal case where a plane wave is incident on the 
v evlinder (the / vector is to be perpendicular to the 
cylinder AXIS . The ‘nits possible to regard M r, 6 
s being proportional to the surface current excited 


. r ™ “ oa oor “- “ ” on the evlinder In other words. t he axial slot can 


FiecuRt ] Implitude of cylinde $ pace facto 


M(x) 


s 
% 
= 


: 


teas 


» 1s 


Azimuth Angie @ be regarded as a receiving antenna and the soures 

taken to be at infinity 
Ficure 2. Phase of cylinder space factor for narrow arial slot, In view of the reciprocal nature of the problem, It 
B. verti« ile is shifted h eu is desirable to compare the present results with those 


256 








he , of Fock [11] for the currents excited on a curved sur- 
he |_| face by an incident plane wave. Fock first considers 
| _}_—-+—J | the diffraction of plane wave by a paraboloid and, in 
particular, focuses his attention on the tangential 
magnetic field in the penumbral region (near the 
boundary of light and shadow). Then, after postu- 
\ Xe lating that the surface currents are only dependent 
' \ , on the local radius of curvature, be obtains a repre- 
= . \ sentation for the current distribution that is only de- 
. { \ , pendent on the parameter //d where / is the distance 
= \ \ from the light-shadow boundary and d is the width 
\ \ — of the penumbral region. Fock then believes that 
\ \ 2 his results are reasonably valid for any shaped bodies 
: \ ° so long as the radius of curvature is always very large 
’ Fs compared to the wavelength. As a crucial test of 
\ Fock’s approximate formula, his results are plotted in 
‘ wo figure 3a along with the evlinder computed data. The 
egreement is very good. There appears to be a slight 
= Y discrepancy for larger positive VY values which is not 
‘d ‘ entirely unexpected because some of Fock’s restric- 
he ' tions are becoming violated. Nevertheless, consid- 
a ' , erable justification is given to the validity of the 
It T 1 extrapolation of the cy linder curves to larger values 
if * \ oll of x. 
l =e 
| | “TS3 3. An Application 
d An interesting application of the foregoing theory 
wo (4-2 is the calculation of the pattern of an end-fire array, 
" of slots on a gently curved surface. For example 
7 consider the array of 2N+1 parallel and axial slots 
_ FiceRE da, Amplitude and phase ef space factor in para- | on a cylindrical surface of radius a indicated in figure 
nal eee: 1. The angle defines the direction of the observer 
re } | 
le 
r . o 
r 
) 
tt tT tr Observer 
(| +A} 
% 2 2Ne! 
Is slot elements 
f AN VIEW 
(| 
° 
( 
th x \. ; *) / 
= / 
' / 
P 
( END VIEW) \ 
, \ 
\ 
AT 
. ‘ \ A 
t Fictre 3b lymplitude o pace clor in parametric forn Y ; ; 
Fieure 4 Schematic diagram of end-fire array of slots on a 
e From ( , , curved conducting surface 
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with respect to the normal at 0, the center of the 
array. The array and the observer are taken to be 
in the principal plane (i. e., the plane of the paper). 
s is the distance from 0, measured along the cylinder 
surface, to P, the location of any one of the elements. 
The field due to the source at P is then proportional 
to 
e*”K(s) A(s) 

where A(s) is the relative strength of the element, » 
is the projection of s on the ray from 0 to the observer 
and A(s) is the correction factor to veometrical optics. 
In terms of the function .VW/, A(s) is given by 


A(s MCN) for X¥ <0, Oa 
x 
~ MNX)e * for X=0, 6b 
where 
\ rT S ka 
(5-*)+5|(3) 
The factor y is given by y=e sin [#—(s/2a)] when e 


is the chord length OP and because 


it follows that 


Ys ( l 3 ) sin ( p = ) 7b 


subject to s a. 
Remembering that there are 2N-+-1 
spaced at intervals As, it is seen that 


elements 


2NXAs=L 
The 


where Z is the length of the array. source is 
now written 


K(s)=e7"* 


where T is the propagation constant of the exciting 
wave along the array. For the present purpose 


3 Lamk, 


where m is a constant real number. 
The total field of the array is now proportional to 


. 
>) e*e-* \(s 
nm —N 

where s=n-As with n N, V+1. 1, 0, 


Introducing dimensionless pa- 


Lo «ao ml, A. 


rameters, defined by 


AS=kAs and S=nAS 


it follows that the pattern 7(#) of the array is given 
by 


T() mS] A(X), (Sa) 


where 


Ss S 
34 (kay?) *" (#3 (hay ) 


*)ti lls] 


When ka tends to infinity such that the surface jg 
effectively flat, the pattern becomes 


and 


= 
ee, 
— 
a 


\ 
1(& T(@ >, kG Qa 
os 
where V(S S sin ® and since S=n-AS, 
on ( on l ef CN TDA 
A - K Ob 


where A=AS[sin 6— 7m]. 

According to the analysis of Hansen and Wood- 
vard [12], the maximum gain in the end-fire direction 
i. e., &€=90°) is obtained when m~(N+2)/N. 
Then, taking the electrical spacing between the 
elements as 45° or AS—7/4. the array has 65 ele- 
ments, V=32 
wavelengths. The pattern 7)(®) in this case is only 
defined for the region 90°< b6< 90°. The main 
lobe and the first side lobe are shown plotted in 
figure 5 normalized such that the maximum field is 
0 db. Using the same values of N and AS, the cor- 
responding pattern J(#) is plotted for ka=200. 
The main lobe is seen to be somewhat changed in 
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Figure 5. Pattern of ¢ nd-fire array of slots on a curved surface. 
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form, particularly in the directions near &=90° 
which is the tangent plane of the center of the array. 
The decibel level of the field in the diffraction zone 
@>90°) is seen to decrease linearly with angle. 
There are no lobes formed in this region. It is of 
interest to know how the curvature of the array 
‘tself would modify the pattern. Therefore, for 
sake of comparison, the pattern 7() of the end- 
of isolated elements on a circular are 
ka=200) is also plotted in figure 5. The working 
formula here is identical to equation (Sa) with 
A(X) being replaced by | because diffraction effects 
are no longer present. The pattern is almost identi- 
eal to that for an end-fire array where elements are 
on a straight line. in this case, are formed 
on both sides of the main beam although they are 
not quite symmetrical about it 

The pattern 7(®) is also computed for ka=100, 
and is shown plotted in figure 6 along with the curve 
for ka and the curved array of isolated elements. 

It would seem that an end-fire array of slots on a 
curved surface has some rather Interesting properties. 
The main effects of the non-flatness of the surface is 
to tilt the maximum of the main beam up slightly, 
and to extend the total width of the main beam into 
the shadow region. There is no evidence of lobes 
on this side of the beam, although at large angles 
6 approaching 180°), there is the possibility that the 
back lobe from the end-fire array could creep around 
the surface in the other direction and form lobes. 
This effect trom a practical standpoint, however, is 
completely insignificant for ka>100. 


fire array 


Lobes, 
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4. Concluding Remarks 


It is seen that the radiation pattern of a flush- 
mounted antenna is influenced, to a considerable 
extent, by the radius of curvature of the surface on 
which it is mounted. In most cases, the main beam 
is tilted upward away from the surface, although the 
total width of the beam can be broadened into the 
shadow. Similar phenomena had been observed in 
the calculated patterns of slots on a conducting half- 
plane [13]. It would appear that any attempt to 
reduce the side lobes of an end-fire array by using a 
diffracting edge or surface will lead to a broadening 
of the main beam. It is also quite apparent that 
the resultant pattern of a flush mounted antenna is 
not simply obtained by multiplying the free space 
pattern by the pattern of a single slot on the curved 
surface. Such a process would produce a pattern 
with lobes in the diffraction or shadow region which 
do not exist in reality for a nonclosed surface. 
Furthermore, this ‘‘multiplication” technique would 
not predict the broadening of the beam in the tangent 
plane. 


The author thanks William Briggs for his assistance 
in the calculation of the patterns in figures 5 and 6. 
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Standards for pH Measurements from 60° to 95° C 


Vincent E. Bower and Roger G. Bates 


The NBS standard pH seale is defined ir 
pH values (pHs 


have heretofore been assigned only in the range 0 


1 terms of several reference solutions whose 
to 60° C. Recent 


measurements of the standard potential of the silver—silver-chloride electrode at temperatures 


above 60 


now permit this scale to be extended to 95 


C, 


Electromotive-force measurements of cells with hydrogen and _ silver—silver-chloride 


electrodes and containing buffer-chloride soluti 
60° to 95° C. 


ons are reported in this paper for the range 


The assignment of pHs values to the following five reference solutions is 


described: 0.05-m potassium tetroxalate; potassium hydrogen tartrate (saturated at 25° C); 
0.05-m potassium hydrogen phthalate; 0.025-m potassium dihydrogen phosphate; 0.025-m 
disodium hydrogen phosphate; and 0.01-m borax. 


1. Introduction 


The Bureau has recommended a standard pH 
scale defined in terms of several reference points. 
The primary standards which fix the points on this 
scale are six in number and cover the pH range 1.6 
to 12.5 at 25° C. Their compositions are as follows: 
4 0.05-m solutionof potassium tetroxalate (pHg, 1.68 
1]'); a solution saturated with potassium hydrogen 
tartrate at C (pHs 3.56 [2]); 0.05-m potassium 
hydrogen phthalate (pH 4.01 [3]); a solution 0.025-m 
with respect. to both potassium dihydrogen phos- 
phate, and disodium hydrogen phosphate (pHg 6.86 


OF J 


$]): 0.0l-m borax (sodium tetraborate, pHs 9.18 
5}); and a solution saturated with calcium hydroxide 
at 25° C (pHs 12.45 16] , 


In earlier published work, the pus values of five 
of these reference standards have been determined 
at intervals of 5 deg from 0° to 60° C, permitting 
the standardization of pH equipment over a rather 
range of temperatures. At the time when 
most of these values were obtained, standards at 
higher temperatures were precluded by a lack of 
data for the standard potential of the silver—silver- 
chloride electrode, or other reference electrode, above 
60°C. A need on the part of commerce and indus- 
try for accurate standards for measuring pH at 
temperatures up to 95° C, however, led the authors 
to undertake a determination of this important 
electrode potential in the range 60° to 95° C [7]. 
The pHg values of all but oneof these standards have 
now been extended to this high range of temperature. 
Calcium hydroxide has a negative temperature 
coefficient of solubility, and a separation of solid 
material was sometimes observed when a solution 
saturated at 25° C was heated to temperatures 
above 60° C. Furthermore, the poor stability of 
those cells in which precipitation did not occur made 
precise measurements extremely difficult. For these 
reasons, pHs values have not been assigned to the 
calcium hydroxide standard in this extended tem- 
perature range 


wide 


1 this paper 


261 


2. Method of Assigning pH, 


The procedure by which the standard pH values 
have been assigned to the standard reference solutions 
has been set forth in earlier publications [1 to 6, 8, 9]. 
There are three steps in the procedure: (1) Measure- 
ment of the emf of the cell Pt; H.(g), standard 
reference solution, KCl (m), AgCl; Ag. This cell 
contains the reference solution in which a soluble 
chloride (usually potassium chloride) is present at 
one of three selected low, accurately known, molali- 
ties (m). From the emf (/), the standard potential 
of the silver-silver-chloride electrode (— E° [7]), and 
natural constants, values of the quantity pwH are 
calculated by the formula 


(E—E°)F 

pw = 53026RT tS ™ ) 
in which F is the faraday, R the gas constant, and 
T the temperature on the Kelvin scale; (2) the effect 
of the added chloride on the value of pwH is removed 
by extrapolation to m=0, in order to obtain the 
limiting value, pwH®; (3) the standard pH values. 
termed pHg, are calculated from pwH® and a con- 
ventional definition of the activity coefficient of 
chloride ion, fe: 


pHg=pwH?® + log fey. (2) 
The single ionic activity coefficient, fo, is related in 
some arbitrary manner to mean activity coefficients, 
which are either experimentally determined or cal- 
culated from equations, like that of Debye and 
Hiickel [10], believed to be valid for dilute solutions. 
For this step it is necessary to know the ionic strength 
of the solution. 


3. Experimental Procedures 


In view of the high vapor pressure of water in the 
cell solutions at these elevated temperatures, each 
cell vessel was provided with a saturator consisting 
of three chambers in series [11], in addition to the 
single saturating chamber that was a part of the 
vessel itself [12]. The hydrogen gas bubbled through 





a total depth of nearly 30 cm of solution before 
entering the hydrogen electrode compartment of the 
cell. By this means, vapor-liquid equilibrium was 
believed to be established at all of the temperatures. 
The smallest volume of hydrogen sufficient to 
separate the solution from the rubber stopper holding 
the electrodes was allowed to remain in the silver 
silver-chloride electrode compartment at the time 
the cell was filled. 

The constant temperature bath was filled with 
a commercial heavy mineral oil of U. S. P. grade. 
Temperatures were measured with a platinum re- 
sistance thermometer, and are believed to have 
been fixed at the nominal values within +0.05° C. 
In order to correct the observed emf to 1 atm 
partial pressure of hydrogen, the vapor pressures 
of the solutions were taken to be the same as those 
of pure water. The pressure corrections at 60° to 
95° C, in millivolts, have been tabulated elsewhere 
[9]. 

The cell solutions were prepared from weighed 
amounts of the buffer materials (NBS Standard 
Samples), and from potassium chloride solutions of 
the desired molality. Dissolved oxygen was _ re- 
moved by bubbling pure hydrogen through the solu- 
tions for 2 hr before the cells were filled. Concentra- 
tions were expressed on the molal scale (moles per 
kilogram of water) exclusively, inasmuch as the mo- 
lality remains unchanged when the temperature of 
the solution is altered. The dilution effect on the pH 
of these standard solutions is small. Hence, in 
practice it is permissible to regard the buffer con- 
centrations given in the tables as either molarities 
or molalities. The error exceeds 0.001 pH unit only 
for the tetroxalate solution, where it may amount 
to 0.002 unit. 


ry’ ° 
The preparation of the electrodes has been 
TABLE 1. 
Reference solution rempera 
ture 
Cc 
4 
70 
0.05-m potassium tetroxalats SU 
90) 
a5 
i) 
70 
Potassium hydrogen tartrate, saturated at 25° C sf) 
90 
QS 
mm 
70 
0.05-m potassium hydrogen phthalate sO) 
“wt 
| ) 
th) 
0.025-m potassium dihydrogen phosphate, 0.025-m diso- 70 
dium hydrogen phosphat« x) 
a“) 
| Ns 
i) 
| 70 
be 
wm) 


0.01-m borax | 
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Number of cells for mxci 


described elsewhere [9]. The hydrogen electrodes 
used in the phthalate solutions were coated with 
palladium black to minimize reduction of the solute 
In other details, the procedure followed closely thay 
of the study of hydrochloric acid solutions at 69° 
to 95° C. [7]. 

The duration of stability of these cells is eop. 
siderably reduced at elevated temperatures, partic. 
ularly when the cell solutions are alkaline. For this 
reason, the cells were brought to initial equilibrium 
at room temperature, and the measurements at 
60° to 95° C were completed within 8 hr after the 
temperature of the bath was raised. 


4. Results 


The measurements were made at three low 
molalities (m), namely 0.015, 0.010, and 0.005, of 
potassium chloride. It has been found that the 
values of pwH for a given buffer solution, computed 
by eq (1), usually vary linearly with m_ below 
m=0.015. Although most of the measurements at 
high temperatures were of insufficient precision to 
provide an independent test of linearity, neverthe- 
less, a straight line, located by the method of least 
squares, was drawn through the values of pwl. 
The intercept of this line 


pwHh= pwH°+am (3 


with the axis at m=0 is pwH®, and the slope of the 
line is designated a. 

The experimental data, which were rather numer 
ous, are summarized in table 1. In the next to 
the last column is given the value of ¢, the standard 
deviation of a single ‘‘observed”’ value of pwH from 
the calculated value, and in the last column is found 
o;, the standard deviation of pwH?®, the intercept. 


Summa? j of data and calculate d values of pwH 


Constants of eq (3 Standard deviations 


0.010 0.015 pwlt a a o 

2 2 1. 827 0. 28 0. 003 0. 003 
2 2 1. 849 35 003 003 
2 2 1. 877 34 O02 002 
2 2 1. 904 16 O12 O11 
2 2 1. 919 ol O11 010 
{ 2 3. 643 42 Oo 002 
1 2 + 664 20 003 002 
4 2 3. 608 “i oo2 002 
t 2 3. 738 25 003 oo2 
t ~ $. 767 1. 56 005 004 
i 5 4.175 05 005 003 
; $ 4.219 03 OO7 005 
; ; 4. 259 ol 006 005 
3 $ 4. 301 + 01 006 005 
1 4.331 04 005 003 
{ 2? 6. 048 o7 005 004 
4 2 6, 962 4 O04 003 
4 2 6.979 2 004 003 
{ 2 7. OO1 + OF 005 003 
‘ 2 7.014 03 004 003 

; 9 026 4 005 003 

; &. G00 2 004 003 

; &. OAS i) 003 002 

3 &. O20 1.19 005 003 

3 &. 809 l.@é 012 Oo9 
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5. Calculation of pH. 


Several different conventions for the estimation of 
the activity coefficient of chloride ion, fe, in eq (2), 
have been used in assigning pHg values to the six 
NBS standards. Inasmuch this coefficient has 
no independent thermodynamic definition, one con- 
vention cannot termed “more correct”? than 
another. In fairness, all of these may be called 
“reasonable”, and it is fortunate that all lead to 
identical pH values at very low concentrations and 
to values that differ by only about +0.01 unit at an 
ionic strength of 0.1 [8, 9]. The assignment of pHg 
values with an accuracy of +0.001 unit must await 
the universal acceptance of a single convention for 
the estimation of the individual ionic activity coeffi- 
evaluation of electrochemical 


as 


be 


cients needed in the 
data. 

Mean activity coefficients of the electrolytes 
present in buffer solutions are ordinarily unknown; 
hence, assumed relationships among the activity 
coefficients of the ions in mixtures are not 
usually a fruitful means of evaluating fc;. It is just 
as satisfactory, and more convenient, to estimate 
fy by the Debve-Hiickel equation 


these 


Ayu (4) 


log fe 
Ba*, bi 


where A and B are constants for the water medium 
at a partic ‘ular temperature, » is the ionic strength, 
and a* is an “‘ion-size parameter’. When the latter 
is given a value of 4 to 6, the right side of eq (4) 
reproduces the logarithm of the mean activity coeffi- 
cient of many uniunivalent strong electrolytes very 
successfully. 

In order to use eq (4) to calculate pHg at elevated 
temperatures, it is necessary to know, within about 
10 percent, the ionic strength of the buffer solution 
at 60° to 95° C. The evaluation was accomplished 
as follows for the five buffer solutions: 

1. 0.05-m potassium tetroxalate: in the earlier 
work [1], the ionic strength : the 0.05-m solution 
was found to be 0.077 at 0° C and 0.074 at 60° C. 
This change was not only sm: all but approximately 


linear; hence, the values of the ionic strength at 
60° to 95° C were obtained by extrapolation, and 
an average value of 0.073 was selected. 

2. Potassium hydrogen tartrate, saturated at 


jr0 


C: the equation for the calculation of the ionic 
strength is given in an earlier paper [2]. The ratio 
of the dissociation constants, K,/K,, needed in the 
calculation, was estimated with the caued ac- 
curacy from the known ratio in the range 0° to 60° C 
[13] by extrapolation to the higher temperatures. 
In this way, it was shown that the ionic strength 
of the tartrate solution saturated at 25° C is prac- 
tically constant at 0.04 in the temperature range 
5° to 95° C. 


3. 0.05-m potassium hydrogen phthalate: the ionic 
strength of this solution between 25° and 60° C was 
calculated from the compositions listed by Hamer, 
Pinching, and Acree [3], and found to change only 
0.3 percent when the temperature is increased by 

> C. ‘Lhe ionic strengths were plotted, and the 
curve extended to 95° C. An average value of 
0.053 was chosen for the range 60° to 95° C. 

4. 0.025-m potassium dihydrogen phosphate, 
0.025-m disodium hydrogen phosphate: inasmuch 
as negligible amounts of free hydrogen and hydroxy! 
ions are present in these solutions, the ionic strength 
is constant at 0.1 from 0° to 95° C. 

5. 0.01-m borax: the ionic strength is 0.02 and is 
unchanged by the increased degree of hydrolysis at 
the elevated temperatures. The formation of poly- 
boric acids is not believed to be sufficiently extensive 
at this low concentration and relatively high pH to 
affect the ionic strength appreciably [5]. 

These values of the ionic strength were used to 
compute pHg by eq (2) with the aid of the conven- 
tional definition of f-; embodied in eq (4). The ion- 
size parameter a* of the latter equation was given 
two values, 4 and 6, and the values of pHg given in 
table 2 are the mean of the two results to the nearest 


0.01 unit. The pHg of the six NBS standards for 
all the temperatures at which an assignment of 
pHs has been made are plotted in figures 1 and 


and are included in table 2 for the convenience of 
the reader. 

The reference solutions should be preserved at 
room temperature, where they are known to be stable 
for several weeks, and portions heated to the ele- 
vated temperatures as needed for the standardization 
of pH equipment. 


pu, of the six reference solutions* from 0° to 95° C 


TABLE 2, 


| Ca(OH)», 


0.05-m KH 0.05-m 0.025-m 
t K tetrox-| tartrate, KH KH2PO,, 0.01-m (saturated 
alate aay: “d phthal- 0.025-m borax | at 25° C) 
at 25° C ate NaoHPO, ¥ 

0 1. 67 4.01 6. 98 9. 46 13. 43 
5 1. 67 4.01 6. 95 9. 39 13. 21 
10 1. 67 4.00 6. 92 9. 33 13. 00 
15 1. 67 4.00 6. 90 9. 27 12. 81 
a 1. 68 4.00 6. 8&8 9. 22 12. 63 
25 1. 68 3. 56 4.01 6. 86 9.18 12. 45 
30 1. 69 3. 55 4.01 6. 85 9.14 12. 30 
35 1.69 3. 55 4.02 6. 84 9.10 12.14 
4) 1.70 3. 54 4.03 6. 84 9. 07 11. 99 
45 1. 70 3. 55 4.04 6. 83 9.04 4 11.84 
0) 1.71 3. 55 1. 06 6. 83 9. 01 11.70 
sa) 1.72 3. 56 4.07 6. 84 8. 99 11. 58 
“oO 1.72 3. 56 4.09 6. 84 &. 06 11. 45 
70 1.74 3. 58 4.12 6. 85 S. OS e 
gO) 177 3. 61 4.16 6. 86 &. 8Y 
a0 1.80 3. 65 4.20 6h. SS 8. 85 
95 1. 81 3. 68 4.23 6. 89 8. 83 


«In the preparation of the buffer solutions, it is usually more convenient to 
dissolve the buffer substance in water and dilute to a given volume in a volumetric 
fiask than to add a known weight of water. As noted above, the pH, of solutions, 
whose molarities at 25° C are the same numerically as the molalities given in 
table 2, will not differ appreciably from the values listed. 
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A Numerical Study of Dedekind'’s Cubic Class Number 


Formula 


Harvey Cohn ! 


In this paper an analytic formula of Dedekind is used to compute class numbers for a 
sample of pure cubic fields, employing the National Bureau of Standards electronic com- 


puter, the SEAC 


ness of combining integral and decimal arithmetic. 
than those of pure cubie fields previously studied by 


much greater magnitude 
hand computer methods 


1. Introduction 


The discovery of fields with high class numbers * 
has become an end in itself {10}... There are gener- 
ally two possible procedures: One is the algebraic 
method, consisting, operationally, of a certain sorting 
of values of polynomials according to factors; the 
second is the analytic method, based on zeta-function 
residues. The second method is of special interest 
when applicable, because of the scarcity of t\ pes of 
field for which a usable formula emerges. 

In 1900, Dedekind applied both methods to the 
“pure” cubic field, or the field formed by adjoining 
\/* to the rationals [4]. The largest class number h 
that Dedekind achieved was h=9 for M=91. In 
1950, Cassels published a table of class numbers |1| 
with the larger A= 12 for \J/=43, using the algebraic 
method. To achieve a substantially larger A it is 
necessary to take larger values of 1/7. Furthermore, 
Dedekind’s analytic formula is the simpler method 
to program for a machine for certain special VW, 
namely, those that equal (s l)/e* for s and e 
integers. A sample of 21 cases were run, vielding 
six class numbers A between 18 and 27, as well as 
other smaller ones (ineluding the cases cited above). 

As an incidental result the quadratic fields (whose 
relation to the pure cubic field is similar to the rela- 
tion of rational numbers to the quadratic field) are 
studied in some numerical detail 

The computation uses a mixture of integral and 
decimal arithmetic like the earlier work with Gorn 
3], but, here, more fortunately, round-off provides 
no difficulty, even far beyond the range covered. 
The work was begun about July 1956 and the runs 
were completed November 1956 


2. Dedekind’s Formula 


The first algebraic investigations of class number 
were made by Markoff [S!, but the first explicit 
“closed”? formula for class number was derived by 


Dedekind. The entire numerical study undertaken 
Present address, Department of Matt Washington University, 
“t. Lou Mo 
A general background in algebraic numb« neluding quadratic forms 
sumed (see [7 
Figu n brack« licate tl t the end p 
L34677 7 


The computation is one of major magnitude, and it iliustrates the useful- 


The class numbers obtained are of 
means of 


here stems from his formula, which we give without 
proof (in an equivalent form): 


ITH (wy) 
ITH (wp) 


h 


(2.1) 


Here h is the much desired class number of the 
real field RCM), and e(>1) is the (often elusive) 
fundamental unit, which is discussed later. The 
function H(z) is defined for z a positive imaginary 
quantity (Imz>0), as the (positive) function 


H(z) =n(z)n(—Z)|[e(e—2)]*, (2.2) 


where n(2) is the (analytic) Dedekind eta-function [6]. 


n(z) =exp (wiz/12) II (1—exp 2mizn), (2.3) 
n=! 
and = denotes the complex conjugate of 2. The 


symbols wy, wr, are rather complicated and are 
treated in the next section, but more immediately 
relevant to the cubic field itself is the following fact: 
If \f is decomposed by positive a, 6, ¢, so that 


M=ab’e’, (2.4) 


(so that a and 6 are free of quadratic factors), then 


R(M'*) is generated equally well by (ab?)* or 
(a*b)**. Hence the number pair (a,6) of relatively 


prime integers (not both 1) uniquely determines a 
field AK, ,= RCM"). For this field the conductor k 
is defined as 


k=ab a®—b? =0 (mod 9) i 
(2.9) 
k=3ab a®*—b? 40 (mod 9). 
As a final definition, 
D 3k? (2.6) 


is called the discriminant. Its significance lies in 
the fact that the algebraic integers, in particular the 
fundamental unit ¢, together with the algebraic 
conjugates «’, €’ (the complex conjugates) satisfy 


D<\(e—e’) (e’ —€’) (e—e’)|?. (2.7) 
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3. Quadratic Field 





time being to mention that the function H(w), ae. 





cording to a well-known property of the Dedekind —_ 
Associated with the cubic field in question is | eta-function, is invariant under a unimodular change ing \ 
the quadratic field generated by adjoining to the | of variables; thus a unimodular integral change of thus 
rationals the imaginary cube root of unity p given | variables of Q(z,y) will not affect H(w) : neither will 
by (say) the positive imaginary root of the change w’ ®. Thus let Q(z,y) be equivalent 
to the reduced form But 
p-+p+1=0. 3.1 imate 
Q*(A,Y)=R.V+SXY+ TY: (3.11 
For this the norm function N is defined by 
Its root may be written Se aes 
N(p+qe)=(p+qe)(p+9e*)=P*—pqt+¢g. (3.2 | token 
‘ A ie. see t 
We finally associate with the field A,, and its ™ R , are & 
conductor & the triples of integers A, B, C}, such 
that so that 
AC=k 3 the ] 
0<S/2<Rs (—D/3)#=] (3.12 . 
1I<B<A 3.4 a a 
because, for all 9%’’ forms corresponding to con- | ie 
(A,N(B , C’p) l. 3.0 ductor k, unit 1 
For every such triple of integers we construct the dD IRT +S 3k. (3.13) 7S 
quadratic form : a {} 
. | We are now guaranteed by eq (3.12) a much more | calcul 
Q(z,y) = Nl Arc+ (B+ Cp)y] 3.6 favorable sealing for the computation of H(e). 
Thus 
and consider any prime p, not dividing 3°, repre- 
sented by ((a,y) for integral z and y. Then{A,B,C . ~ - 
is called a residue (or nonresidue) triple or ((z,y) is log H(w rk3 12k) Dlog 1— 2x" cos Sn/R+1n* ' Th 
called a residue (or nonresidue) form according t k=exp(— 2rk3*/R), ee 
whether or not the equation (3.14) | 44 bj 
ange 
x =a*b(mod p 3.4 and «Sexp(— 73"), so that to the degree of round-off ae 
; ; se . used here, it will develop that no more than five | ¢hird 
is solvable for an integer z Then for each triple | terms of the infinite sum (or product) are necessary, nd 
satisfying eq (3.3), (3.4), or (3.5) a form ((7,y) 18 | when the reduced form is used for Q(x.y). i . . 
defined with “‘root”’ on 
’ quoti 
B+Cp 4. Units and Choice of Field | divisi 
—— { oS : i quoti 
5 We henceforth confine ourselves to fields generated 8 
The root » Is an wz or an wy according as the residue by M where, for a positive or negative integer 8, 
or nonresidue character of the form. The product Vl , 41) 
symbol in eq (2.1) should then be self-explanatory. : ; ' = 
It can be shown there are 3k” triples {A,B,C} that T i alia ~s 
generate an wp, and 6k” that generate an wy (not all ie advantage 1s that a very obvious unit Is 4. 
of the wp or wy being distinct), where 94”, the total vA , vi Vv ‘ 2 
number of triples, is given by ” al ie sl ia ;' , (4.2) — 
). 
Ok” —kII(1 3 p)/p). »g where M refers to the real Foot. A fundamental unime 
pik unit is clearly of the type y'/", and the value of m fd 
can be ascertained in the sample treated here, by equiva 
Here 3/ Pp) is the usual quadratic residue symbol, means of well-known results of Dedekind supple- 3.12) 
mented by the lemma (compare [9]): of Y: 
: Lemma. In the notation of (2.5), if k>10s, then | for wl 
| for p=1(mod 3) 
* pl /2 2) 
3/p) 1 for Pp l(mod 3 3.10 a ae 
0 for p=0(mod 3). Proor. Let e=7'/" be the fundamental unit. Then | lf th 
‘ e, as an algebraic number, has a (root) discriminant Euclic 
not less than the discriminant D of K,.». Hence modul 
The properties of the quadratic form Q(z,y) will — 
* See Ni 


It might suffice for the 


7 -—7'9o' >  —sF 
Ds \e—e’ |"e— |" le il a 


dee Se 


eventually come into play. 
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where e’, €’ are the complex conjugate roots belong- 
ing to the same equation ase. But (e’|=|é’|=e' R 


thus 


te°(1+¢7*”)4, (4.4) 


1/2)4 1/2)2< 
Des le+e )*(2e 


But finally, D=3k?, e=n'’", whereas y is approx- 


imately 3s? from eq (4.2). Hence 
abe < 1(38- hes \ | , Ui] - . 1 38°) . (4.5) 
Taking the maximum on the right when m23, we 


see that the corresponding factors in parentheses 
are as follows: 
dk? S 4(38°)-2*(1 4-2! +-2°*) /3, (4.6) 
the last factor being the one where m=3, s ), the 
inequality (4.6) is contradicted when k>10s. Q.E.D 
The list of 21 cases is given in table 1 in order of 
increasing k. Those cases in which the Ha! sensor 
unit was known to Dedekind or Cassels are denoted * 


hy a “D” or “C”’ respectively. In the other cases, 
easily, k >10s. The possibility that »=é is excluded 
by the fact that in these cases the class number. 


calculated by assuming n= e, will turn out to be odd. 


5. Subroutines 


subroutines were 
One should bear 


The following number-theoretic 
programed for this computation: 


in mind that the SEAC is a binary machine with a 
44 bit word and sign bit. The numbers treated 
range arithmetically from 1+2-* to 4—2-*, 1. e., 


the decimal point is fixed between the second and 
third high-order bit. Integers are held in the low 
end of the register (or scaled by 2~*). 

Division of two untegers to produce integral 
quotient and integral remainder. 


(The built-in 
division operation produced only a “decimal” 
quotient. ) 


2. The euclidean algorithm, or solution in ra ig of 
s\¥+tY=(s,t). (5.1) 
x § Reduction ot modules. 5 
t. Multiplication of modules 
5. Generation of coefficients of (2,y) from { A,B,¢ 
of eq |. .6)). 
6. Transformation of coe fic rents of ()(r,y) under 


unimodular change of variables. 

Reduction ot (a YY) to the properly or improperly 
Pita Hor van reduced form ()* ry) subject to conditions 
3.12). To do this we op find the minimum value 
of Q(2,y) by scanning all of the finite set of points 


for which O(r,y)<Q(U1,0) y _ or for which 

[Ar+ (B+-C/2)y}?+3C*7?/4< A’. ».2) 
if the minimum value is ((z»,y%)=R, then the 
Euclidean algorithm generates an integral uni- 
modular transformation of coordinates, (or finds 


* See headnote to table 1 
' See section 7 


12,42), 80 that dy%—ayyo=1. Thus in the new 
coordinates, 
oe, & | Ay Y’ & 
(5.3) 
YoX'’ +dnY'’=y 
V(r,y) becomes 
QO(X,Y)=RX"+S'X Y’+T’yY”. (5.4) 


Now to minimize the coefficient of X’ Y’, all we need 
do is perform the transformation 


X’=X+9Y 
(5.5) 
"=f, 
choosing g [S’/2R], (the integral part), and we 


obtain Y*(X,Y) of eq (3.11) and (3.12). 

8. Determination of the residue character of a quad- 
ratic form. This subroutine finds a small® prime p 
represented by Q*(X,Y) and not dividing D. (The 
procedure is similar to the search for the minimum 
of V(a,y) by eq (5.2).) Then the eq (3.7) is tested 
with +2 and —z, where z runs from 1 to (p—1)/2. 
If x passes through all these values without eq (3.7) 
being satisfied, then, and only then, Q is a non- 
residue form. 

9. The built-in “‘base zero’’ operations. These are 
standard fixed-point subroutines for input con- 
version to binary form and output conversion to 
decimal form, for logarithms, exponentials, sines, 
and cosines, among others. 


Class Numbers 


The subroutines listed earlier, taken together 
with the formulas and concomitant explanations, 
should indicate the basic program for computing h. 
It must be added that (binary) nonintegral quanti- 
ties used, e. in logarithmic computations were 
scaled by 27", and the final decimal values for 
h logio « were expressed to six decimal places. 

The basic procedure is to generate the 9k’’ triples 

A,B,C} in lexicographic order. For each triple we 
compute the corresponding forms Q and (reduced) 
(*, and finally we form (and accumulate) the con- 
tributions (3.14), resulting in the following data (in 
decimal form): 


6. First Run. 


(rr 
B*) 
14 


Input: “Case number’, a, 6, k, D 


Output: For a nonresidue form, eq (3.11), 


R,S 


6 No effort was made to secure the smallest p, although this seems to have 
happened accidentally in the few cases spot checked. The smallest prime would 
have required much more programing, although the effect might have proved 
rewarding 
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For a residue form, 
r,p,R,S 
Final output, 
h log ce, h logy « 
For simplicity, 2 and S, which are never negative, 
were packed into one word, likewise xz (with sign 


removed) and p. To identify a residue form, the 
x, p combination was printed first, with a minus 


sign. The residue and nonresidue forms are omitted 
from the final tabulation to save space. The 
division h=A logye/logye was performed by hand, 
resulting in integral A to five significant figures (see 
table 1). The data from tables of Dedekind and 


Cassels are denoted by “D” and “C’’, respectively. 

During the actual runs, the program was read in 
on wire, and the parameters were read in on teletype 
tape in trios of ‘“‘case number’, a, 6. (Here the case 
number was packed with identifying information to 
check the tape position.) The output was on 
teletype printer. 

The program exhausted almost 900 of the 1,024 
electrostatic Memory spaces. 

The running time varied from 3 minutes for case | 
to 14 hours for case XXI. A subroutine permitted 
the current {A,B,C} and cumulative sum of +log 
H(w) to be dumped every half hour or so for rerun 
purposes. Almost no time went into input-output, 
by comparison. The vast bulk of the time went 
toward the testing of the residue characters, 1. e., 
toward the solution or nonsolution of eq (3.7 The 
problem was brought to a close at case XXI largely 
because time was beginning to become prohibitive. 


TABLI 


“D” and “C”"’ denote units or ¢ numbers apy 


Case Ss V a h ( log 7 é 
I l 2 2 l 1 0. 585159 n 

I] 2 9 l 3. 1 1.096455) 
III 3 28 7 2 1 1.436650 » 
1\ 18 58351 17 l 7 2. 987641 n 

V S 913 19 1 | 3) 2. 283584 ” 

VI 2 7 7 | 1 1. 060137 n 
Vil 3 26 26 l 1 1. 425922 ” 
VIII 19 6860 5 2 7 3. 034650 n 
IX 10 999 7 l 3 2. 476964 » 

X | 63 7 } 1 1. 678966 ” 

XI 4 728 9] 1 | 2 2.385408 | 7» 
XII 23 12168 | 9 | 2 3. 200589 ” 
XIII 5 126 14 3 | 1) 1. 876216 ” 
NIV 7 344 I} | 2 2. 167739 ” 
X\ 5) 124 51 2 1 1. 873900 n 
XVI 1 65 65 | 1 1. 683490 ” 
XVII 6 215 215 | 1 2. 032753 ” 
XVIII 6 217 | 217 | 1 2. 034093 o 
XIX 7 342 38 } | | 2.166895 » 
XX 15 3376 | 422 1 | 2) 2. 829347 n 
XXI 17 914 | 182 1 | 3 2. 938049 ” 


The SEAC has no overflow stops at all, not eyep 
for 0/0, and no automatic error checks. As a substi. 
tute for both, a test was made to see if the discrimi. 
nant of each quadratic form checks with —3} 
In the range covered, no overflow would haye 
occurred, barring machine errors. 


7. Group Structure of Modules 


A good deal of insight into the machinery of 
Dedekind’s formula gained by considering the 
composition structure of the group of quadrati, 
forms in the “order” of discriminant 3k. One 
again, rather than summarize a_ well-known and 


Is 


lengthy theory, we shall consider the modules onh 


in relation to the present computation 
The notation 


M=—[A, +B, ..., Act 


denotes the aggregate of integers in R(p 


+ Bip), 


En, (A 
Such a set is called a module 
B+ Cp) 


section 3 


where n, are integers. 
Every module has a reduced basis like LA, 
identified with the triple (A,B,C) of 
Here A taken the minimum positive rational 
integer in M, and B+ Cp is defined by the conditions 
that €’ the coefficient of  p, 
whereas B is the minimum positive integer na module 
element in which the coe fheve nt of p is C. The enumer- 
ation of triples (.4,2,C'} enumerates modules, of 
course, whereas Q(r,y) is the norm of an arbitrary 
element of the module. Clearly 


iN as 


Is minimum positive 


| 
n Dedekir nd ¢ vor} espectiv 
} ke?’ h log t h CGren- D 
erator 
6 l 0. SS5150 | 1) 2 LOS 
if) | 1. 096455 | D> 3 243 
) 14 2 2. 154975 3 ) 28 588 
17 2 2. YS7T641 | 17 R67 
1) 19 2 3. 425376 3 DD 19 1083 
21 2 ». ISO0412 3 DD 7 1323 
26 } 1. 277770 3 © 26 2028 
30 6 t 55S1LO74 3 DD 20 2700 
7 | 7. 430929 3 Cc 37 $107 
63 6 10. O73793 6 D 63 11907 
4] S 21. 468668 4 8) | 24843 
I) 117 12 19. 203533 6 " 39 11067 
1 2¢ Is 16. 885942 y 126 17628 
1) 129 14 26. 012867 12 S 3 $9923 
IS6 30 16. 865096 y 124 LO3788 
DD 195 24 30. 302812 IS 65 114075 
215 | 28 12. 687799 | 21 215! 138675 
217 | 20 4. 920504 | 27 217! 141267 
342 54. 58. 506160 | 27 342 350892 
122 70 59. 416280 21 $22 534252 
546 | 72 79. 327311 27 182! 894548 
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FH toto -1 or 


pM pl A, B+ Cp|=—|pA, C+ (B—C)p] 
is another module leading to the same GV(r,y), as | 
N(p)=1, 80 is pW. Thus the triples {A,B,C} or | 
the reduced modules [A, B+Cp], with AC=k, 


come in “‘threes,’’ leading to the same form Q(z,y). 

Dedekind considered the multiplicative group | 
formed by the modules under multiplication of basis | 
elements: 


A, B, T C1p\| Az, B, T C’,p)| 
(A,A2, AyB,+ Ai Cyp, AgB, + AzCip, (B\B,—C,C,) 4 
(B,C,+- B,C, —C,C2) p). 


(7.2) 
It can be shown that the product module can be 
reduced to the form [A;, B,+C,p| where A,C,=k, 
provided A,C,= ALC k). This process is closely 
analogous to the composition of forms, or, of course 
the multiplication of ideals. 

For the purpose of performing this multiplication 
of modules we need two subroutines mentioned 
earlier (numbers 3 and 4 of section 5). 

3. Reduction of modules, (from [A), 
B,+Cyp| to the reduced form |A, B+ (Cp]). 
4. Multiplication of modules. This is the repeated 
application of the earlier subroutine for the reduction 

of product (7.2) to [A, B+ Cp]. 

The notation “‘residue”’ (or ‘‘nonresidue’’) module 
ean now be introduced to correspond to ‘“‘residue’”’ 
or “nonresidue’’) form. Dedekind has shown that 
the set of residue modules forms a subgroup of index 
three among the set of all modules. Thus the non- | 
residue modules can be partitioned into co-sets so 
that all modules are of type 


B, + Cip, 


where Ni is a variable residue and MN, is a fixed non- 
residue module. In the notation H(M) for 
H(B+(C'p/A), we can write in accordance with (2.1) 


. TTA (RN, 
ITH (9) 


Now if we take precautions that once R appears, pR 
and pk shall not appear, we have only k’’ factors in 
numerator or denominator. 

Now (7.3) is very close to Dedekind’s original 
formula. (Dedekind used quadratic forms and not 
modules to generate the products in (7.3).) 


IS 


8. Ambiguous Forms and Modules 


lo appreciate the relation between modules and 
forms we must realize that not only do M, pM, p?M 
correspond to the same form but 9 and WM can 
correspond to the same form. Here MM, the conju- 
gate module is defined by 


M=[A, B+Cp]} } 


M—=[A, B+ Cp’. 


The difficulty really arises when M=M, (or pM=M, 
etc.). (The second possibility is really no different 
as it implies p*3 is then its own conjugate.) These 
are the ambiguous forms of classical number theory. 
Generally MM=|[1, kp] the unit element, or principal 
ideal so that if M—WM, MY is a unit element and M 
must belong to any subgroup of index three. In 
other words the self-conjugate modules are all 
residue modules. 

To give an ample illustration let us take case VI, 
s +, k=21, and list the 9k’’=18 modules in lexi- 
cographic order, as does the computing machine, but 
with the interrelations in question clearly indicated. 


List OF MODULES FOR k=2] 
(1, 21p]=M, (7, 5+3p]=M, [21, 7+ p]—pM, 
(3, To] =M, [7,6+3p]=M, [21, 9+p]—pM, 
[3, 14+7p]=M, (21, p]l—=pM, (21, 13+ p]l=p’°M; 
[7, 30] =pM, [21, 1+p]=p?M, (21, 15+p]=p7M, 
[7, 3+-3p]=e'M, [21, 4+p]=e'M, [21, 16+ p]=p'M, 
[7,44+3p]=M,  [21,6+p]—pM, [21, 18+p]—pMbo. 


These modules form a group with identity element 
Wy. Rather than construct the whole group table 
we might identify M,,eM,,e?M, as one element, and 
it is easy to see that the group reduces to an abelian 
group of order 6 given by 


My, Mr, Ms, D—=— DME, M.— MMs, DE— DEM, (8.2) 
with the generators YU,, Mi, satisfying 
M,= DE= MF. (8.3) 


We see Ny, Ms; are the residue subgroup, for clearly, 
Mio= Meo, M,= Wes. 
The corresponding quadratic forms, (grouped after 


Dedekind with their “negatives’’), are ‘ 


Mo: 2try+331yY=—Q(z,y) ) 
M,: Pa? + 3ry4+-37y7=Q,(z,y) 
_ (8.4) 
Wi: 132° + 9ry+ 27y°=Q2(z,y) 
Ws: 192" + lleyt 19y°=Q3(z,y). J 





Forms for Pi, and WM, are ambiguous (or self-conju- 


gate). Note that ab?=7,so0 that for the primes given 
331=@,(0,1), 10° =7(mod 331) 
19=Q,(1,0), 4°=7(mod 19) 


The negative of O(7,y) is found by replacing S by — Sor, effectively, by keeping 


R and D fixed and replacing S by 2R-—S 
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in accordance with eq (3.7). In the “general” situa- 
tion, of course, the residue modules will not all be 
ambiguous. 


9. Second Run. Residue Group 


A second run was undertaken to simulate more 
closely Dedekind’s original method, which consisted 
of using the structure of the residue group. The 
advantage in so doing is clear if we think of the fact 
that most of the computing time goes toward testing 
eq (3.7) for a solution. The amount of work can be 
cut to about a sixth * if we agree to keep the quad- 
ratic forms 1n storage so as to never test the same 
form twice. Because the memory had been almost 
exhausted by the program, magnetic tape was used 
for storage of forms. 

The second run proceeded as follows. The ma- 
chine generated the modules, testing each in turn 
until it came to (say) 9), the first nonresidue module, 
which it put into memory. (Only the internal 
memory had been used thus far.) Once the machine 
reached this point it returned to the beginning of the 
list of modules and did two things “simultaneously”: 
First of all, it tested each module whose quadratic 
form was not on the tape and loaded on the tape ° 
every form that had been tested. Second, the ma- 
chine took the residue module sav) ® that had not 
been previously tested (and loaded on tape) and 
formed IR loading both quadratic forms on tape, 
while it accumulated log H(2, 2) —log H(M). 

The second run led to an interesting array of 
quadratic forms which we see in table 2. Here R,S 
refer to coefficients of the residue quadratic form 
and R’,S’ refer to the corresponding nonresidue 
quadratic form which is formed by multiplying by 
the module ¥, as explained earlier. 

Unfortunately, here overflow was a serious prob- 
lem in the multiplication of modules, and only cases 
| through XVI, excluding XIV, went to completion. 
The forms for which S= RP or zero obviously correspond 
to self-conjugate modules."° The other self-conju- 
gate modules are marked by ‘‘*’’, and can be spotted 
by the fact that the residue quadratic form and its 
negative do not both appear in table 2 (see footnote 
4) 


§ Unfortunately tape motion produced too much delay is much as 1 min per 
search to permit these benefits to be fu 1 

Here, whenever =I, pM, or pW, one should load both the QO for M and the 
negative Q (for Pi) on the tape to prevent both Q and its negative from appearing 
is residue forms This fact had been erlooked and had to be corrected by in 
indesirable manu ntervention in the ! 


See footnote to tabl 
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10. Concluding Prospects 


The computation consummated here represents 
one phase of what really amounts to a larger task. 
The computation of units should also be programed, 
possibly after Voronoi’s algorithm [5], so that A 
can be found even when M does not have the for- 
tuitous form s*+-1. The storage of the previously 
tested forms in a larger high-speed memory would, 
of course, have speeded up computations, but the 
testing of congruence (3.7) leaves much to be desired 
in the direction of efficiency. To carry the groupings 
of modules and forms on to larger k, a double pre- 
cision multiplication of modules would have to be 
provided. Lastly, it would be desirable to find the 
class structure algebraically [2], which is not alto- 
gether beyond the capacity of the SEAC. 

The class structures, as well as the group struc- 
tures, of the modules shall be the subject of further 
theoretical investigations. 


The author thanks Olga Taussky for her aid and 
encouragement. 
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Impact Properties of Slack-Quenched Alloy Steels 


Melvin R. Meyerson and Samuel J. Rosenberg 


A method was developed for slack quenching impact specimens so as to produce uniform 
and accurately controlled hardnesses and microstructures at any location in the cross sec- 


tion. 


Steels of the 8Ixx type, without and with boron, and varying in hardenability, were 


studied to evaluate the effect of slack quenching upon the impact properties. 

Slack quenching was shown to be detrimental to the impact properties of the steels. 
The injurious effects were found to be dependent on the carbon and alloy content of the 
steels, the degree of slack quenching, and the amount of tempering. 

When any two steels were slack quenched to the same hardness, the higher carbon 
lower alloy steel had impact properties inferior to those of the lower carbon higher alloy steel. 
Tempering usually improved the properties of both and reduced the differences between 


them 


The higher the hardness to which any steel was slack quenched, the lower were its impact 


properties 


When tempered to the same hardness 


however, the steel originally slack 


quenched to the higher hardness had the better impact properties. 


1. Introduction 


Despite the fact that slack-quenched structures are 
frequently found in heat-treated steels, quantitative 
evaluation of the effect of slack quenching upon the 
mechanical properties of such steels has never been 
adequately defined. This is due, not to lack of recog- 
nition that such data are needed, but rather to the 
lack of heat-treating techniques that can accurately 
control the desired slack-quenched microstructures. 

In an effort to overcome this deficiency, many 
researches have evaluated the effect of nonmarten- 
sitic decomposition products of austenite upon the 
mechanical properties of various steels. With but 
one exception, the amount and nature of the non- 
martensitic products have been controlled by the 
use Of isothermal transformation techniques. Al- 
though much excellent work of this nature has been 
earried out, the results obtained are not necessarily 
indicative of the properties of slack-quenched steels 
in which the amount and nature of the nonmarten- 
sitic products may vary considerably, depending 
upon size effects and the rate of cooling. The re- 
ported exception is the work of Vajda and Busby 
1,2], who end quenched 7-in. diameter rounds and 
tested specimens cut from sections parallel to the 
quenched ends. Although the technique developed 
was new and useful, it suffers from two disadvantages 
in that the mass of steel needed for development of 
the slack-quenched specimens is considerable and 
the cross sections of the test specimens so obtained 
cannot be uniform because finite differences in cooling 
rates must exist. These authors attempted to mini- 
mize the latter by using subsize specimens, but the 
fact remains that finite changes in cooling rate, with 
corresponding differences in microstructure, exist 
over the cross sections of the specimens so obtained. 

The project described in this paper is concerned 
with the effect of slack quenching upon the impact 
properties of steels, data which were desired by 
Watertown Arsenal. The problem posed was: 
Given a pair of selected steels whose Jominy (end 


—_.. 


Figures in brackets indicate the literature references at the end of this paper. | Of 8 was assumed for all steels. 


quench) hardenability curves crossed each other, the 
intersection of these curves represent equal hardness, 
obtained at equal cooling rates. The resulting mi- 
crostructures, however, could be expected to be 
dissimilar. Information was desired as to the effect 
such slack quenching would have upon the impact 
properties of the steels involved. 

The present paper describes a simple technique 
for accurately controlling the hardness and micro- 
structure of slack-quenched impact specimens, and 
gives the results of the effect of slack quenching, 
both with and without tempering, upon the impact 
properties of a series of 81xx steels, both without and 
with boron. Specifically, a comparison is made 
between the structure and impact properties of 
various pairs of steels, slack quenched to the hardness 
corresponding to a point of intersection of their 
Jominy curves (the steels were intentionally designed 
to give this intersection). A similar comparison is 
made of the structures and impact properties of 
pairs of the same steels treated with boron and slack 
quenched to obtain a hardness equal to that at the 
point of intersection of their Jominy curves. A 
third comparison is made between pairs of the boron- 
free and boron steels. A few supplementary tests 
were made on some of the same steels as fully har- 
dened and then tempered to the same hardness 
levels as the slack-quenched steels, or the slack- 
quenched and tempered steels. 


2. Experimental Steels 


2.1. Design, Melting, Fabrication, and Analyses 


Using 8140 as a base, it was planned to study the 
effects of slack quenching upon similar alloy steel of 
higher and lower carbon, with the alloy content so 
modified that the Jominy hardenability curves of 
the various steels intersect. Jominy hardenabilities 
were calculated in accordance with the procedure 
described in the Metals Handbook [3], using the 
values and multiplying factors given therein. For 
purposes of calculation, an ASTM grain size number 
The validity of this 
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TABLE 1. Chemical compositions and hardening temperatures of the steels 
Percent by weight Hardening 

Steel No tempera 

ture 

Cc Mn P Ss Si N Cr Me \ B 
I 
+ 4 0. 49 0. 87 ( 0. 22 0. 30 0. OS 0.05 Py 1 ww 
BR-21 | . . f None \ 
BR 2» f st) ‘1 0. 005 0. 029 IS 24 ; ( On 1 0.0017 1 Ow 
BR-17_..)) ” 7 . = . » f None \ 
BR-IS8 f 44 87 005 026 27 $2 it) lt 0 1 0.0019 i l 0) 
BR-1 | . { Nome ) _ 
BR-2....\J 5 " ” ” = a _ 1 0.0022 = 
BR-25 | _ . . f Nome } 
BR-26 j 3 io 2s 2 6S l vo 1 0.0021 j 1, C00 
BR-19 | . ‘ - . . : f None \ 
BR ” j 2s ”) OO 027 i) 5 HZ 21 OS f 0.0017 1. 600 
, 

+. ~ ; 19 ” OO: 24 $5 63 St) $2 OS a 1, #2 
BR-29 | on ’ f Nome \ - 
BR-30 { 19 ) $3 OS St) 19 O06 1 0.0024 1, 62 


assumption was borne out by subsequent metallo- 


graphic examinations. 


Beginning with the calculated hardenability of 


8140, using the composition corresponding to the 
middle of the range (0.41% C, 0.85% Mn, 0.28% Si, 
0.30% Ni, 0.438% Cr, 0.11% Mo), hardenability 


curves were calculated for steels at carbon levels of 


0.20 percent, 0.30 percent, and 0.50 percent with 
manganese and silicon constant at 
8140 and with nickel, chromium, and molybdenum 
increased for the two lower carbon levels and de- 
creased for the higher carbon steel. In changing the 
alloy content, the Ni-Cr-Mo ratio was maintained 
constant at the same ratio as exists in 8140. Several 
curves were calculated for each carbon level, using 
different totals of alloy content, and from these 
curves four steels varying in carbon content were 
selected for melting. 

All steels were melted in an induction furnace, 
using a magnesia crucible. Each melt was 150 lb, 
deoxidized with aluminum (2 |b/ton) and split two 
ways during pouring with boron (as ferroboron 
added in the ladle to one-half of each melt. The 
steels were poured into hot-topped split molds 13! 
in. high, 2% in. square at the bottom and 3% in 
square at the top. A cast Jominy specimen also 
was poured from each split heat. The ingots, con- 
taining about 45 lb of usable metal, were either 
surface ground or cleaned by tumbling and then hot 
rolled to % in. plate. A 4 in. length (for Jominy 
specimens) was sheared off when the plate thickness 
reached 14%in. All plates were normalized at 1,650° F 
before being machined into specimens. 

The compositions of the steels are listed in table 1. 
Chemical analyses were made only on samples taken 
from the boron-treated halves of the split heats. 

A metallographic examination was made of each 
boron-free steel shown in table 1. Photomicrograplhis 
of these (except BR-29) in the hot-rolled and nor- 


the midlevel of 


malized condition are shown in figure 1. Ferrite 
banding was absent or very slight, BR—1 being the 
worst in this respect, In three of the steels (BR-17, 
19, and 23), the carbides were partially spheroidized, 


2.2. Hardenability 


At least one Jominy hardenability test was made 
on each steel, both as cast and normalized, and as 
hot rolled and normalized. 


The proper hardening | 


temperature, based on carbon content, was used for | 


each steel in making the Jominy tests, and in all 
subsequent hardening. These temperatures 
shown in the last column of table 1. The addition 
of boron resulted in improved hardenability im each 
case, and the amount of improvement varied with 
the composition of the steel (figs. 2 and 3). Gener- 
ally, those steels that had the lowest carbon content 
were improved the most by the addition, an exception 
being BR—1, which was low in manganese and silicon 
The results obtained on both the cast and wrought 
specimens of the same steel agreed closely. How- 
ever, the curves did not, in general, coincide with the 
computed curves; the experimental curves indicated 
lower hardenability than was computed (fig. 4). 

The ASTM procedure for making Jominy harden- 
ability tests [4] specifies that the hardenability test 
bar be quenched from the proper hardening temper@- 
ture for that steel. The graph paper supplied by 
ASTM. for plotting the results of hardenability tests 
indicates that a specific cooling rate (at 1,300° F) ts 
obtained at each increment of distance from. the 
quenched end. Actually this cooling rate can be 
influenced somewhat by the hardening temperature 
However, because the range of hardening temper@- 
tures used in the present study was relatively small 
(465° F), it is believed that variations in cooling 
rate, due to the different quenching temperatures 
used, were insignificant. 
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Fiaure 1. Structures of boron-free steels. 


Hot rolled to 5/8 in, plate and then normalized (1 hr at 1,650° F, air cooled). Etched in 1% nital. X 100. 
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Figure 4 Comparison of computed and actual hardenabilities 


of the steels 


of some 


In order to determine accurately the points of equal 
hardness at equal cooling rates between any two or 
more steels, all of the hardenability curves were 
plotted on an expanded scale as a single master 
graph, and the intersections of the various curves 
determined. For practical purposes, and because 
of the limits of accuracy in hardness testing, curves 
approaching each other within % R, were con- 
sidered as intersecting even if no actual intersection 
occurred. Only intersections at a hardness level 
of 30 R, or higher were used. Four such hardep- 
ability curves, taken from the master graph, are 
shown in figure 5. 

Based on carbon content and provided that the 
hardenability curves cross each other only once, 
there are six possible intersections of two curves 
each, for the boron-free steels, and six for the boron 
steels. There are only six possible intersections of 
two curves each for the boron-free steels versus the 
boron steels, because the presence of boron always 
improved the hardenability. It was found that 
none of the curves of the higher carbon steels con- 
taining boron intersected those of the lower carbon 
steels without boron. For convenience, each com- 
bination will be referred to as follows: 


‘and 0.41°, 
‘and 0.30°; 


Code A—0.50°;, ¢ 
( 
(and 0.20°; 
( 
( 
( 


Code B—0.41% 
Code C——0.30°, 
(‘ode D—-0.50°, 
Code E—0.41°, 
Code F—0.50°, 


and 0.30°, 
and 0.20°, 
and 0.20°, 


— . . o . . 
Although 18 combinations of one intersection were 
theoretically possible, only 15 could be obtained, as 
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Ficure 5. Jominy hardenability curves for four of the steels 


Intersections of the curves indicate points of equal hardness at equal cooling 


rates 
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indicated in table 2. In order to obtain this many | Tasve 2. Selected intersections (+4 R,) of Jominy harden- 


r |}  eombinations, it was necessary to utilize steel BR-1, ability curves 
| which had a low manganese and silicon content. 
- » steels were ex: ine > yr 1° , | Combi- Hardness at) Distance 
I All of the ste . Is were examine d ime tallographic ally nation Steel numbers and percent carbon intersection, from 
is after being slack quenched to the proper hardness code content Re | quenched 
e | jevel at which they were to be tested, and the micro- end, 1/16 in. 
»e +7 ’ "OQ mC yf , ar SS ‘Ves » j 
tS) struc tures at ea h . f the h urdne — levels not d in Boron-free steel versus boron-free steel 
\- table 2 are shown in figures 6, 7, and 8. The struc- 
n | tures varied according to carbon and alloy content, ‘ BR-27 (0.49C) and BR-17 (0.44) 46 5 
} ‘ an 13a wa 1 » aetaale are B BR-1 (0.38C) and BR-25 (0.30C) 49.5 3 
| and the hardness levels to which the steels wer : Swe aa né Bie nee = : 
)- slack quenched, and the two steels of any pair D BR-27 (0.49C) and BR-25 (0.30C) 33 7 
. . 4 ot obtainec 
e quenched to the same hardness usually varied in ; Not obtained 
structure, often markedly. Steels quenched to the 
e higher levels of hardness contained appreciable Boron-free steel versus boron steel 
p amounts of martensite, some bainite, and oceasion- 
; ; ai P . ; —melate « : th <2 AVA \ BR-27 (0.49C) and BR-18 (0.44C) 54 4 
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FIGURE 6 Vicrostructures of pairs of steels, boron-free, slack quenched to hardness levels selected for testing 
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3. Development of Heat-Treating Techniques 
3.1. Controlled Slack Quenching 


A Jominy specimen contains an infinite number of 
planes parallel to the quenched end, each of the- 
oretically constant cooling characteristics (assuming 
no surface cooling), and decreasing in rate of cooling 
- with increasing distance from the quenched end. If 
heat loss through radiation is low in comparison with 
the rate of abstraction of heat from the quenched 
end, impact specimens should be amenable to end 
quenching and so develop an infinite number of 


eren numbered 





(48R,) BR-I7 


BR-2I 
0.46%C 


BR-24 (43 R,) 


and bhoron-free odd numbered). slacl quent hed tothe 


selected for testing. 


l nital x 500 


planes, each uniformly cooled. With this thesis in 
mind, experiments were conducted to evaluate the 
effect of various means of end quenching impact 
specimens so as to develop controlled hardnesses at 
selected distances from the quenched end. 

The root of the V-notch in a Charpy specimen is 
1.08 in. (approximately '%{. in.) from each end. In 
steels of moderate hardenability, the hardness value 
obtained at '%j¢ in. from the quenched end of a Jominy 
specimen usually is appreciably lower than_ the 
maximum obtainable. It was obvious, therefore, 
that any method that was developed for end quench- 
ing impact specimens would have to give markedly 
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FIGURE 8 Vicrostructures of pairs 


Etched with 1° 


faster cooling rates at a distance of '%,_ in. from the 
quenched end than that attained in the conventional 
Jominy specimen, if impact specimens were to be 
slack quenched to hardnesses only slightly below 
the maximum possible. Several methods to speed 
the rate of cooling, involving the use of various 
specially designed shapes and adapters, were tried, 
but none was completely satisfactory. The method 
finally adopted had the advantages of simplicity, 
reproducibility, and economy of material and machin- 
ing. The details of this method are as follows: 

_ The impact specimen is rough machined to approx- 
imately 0.42 in. square to allow for scaling and de- 
carburization and to a length of 2%. in. so that 
there will be some latitude for locating the notch. 


’ : Aho 
(33R,) BR-28 
0.49%C 


of steels, boron-treated, slack quenched to the hardness levels selected for testing. 


nital x SOO 


One end is drilled and tapped to accommodate about 
a \-in. length of a 4-20 screw, 2 in. long. A 
washer, locked in a selected position by two nuts, 
serves to support the specimen in the quenching 
fixture so that a predetermined portion of the 
length of the specimen is immersed vertically into 
a still-quenching medium, which may be brine, 
water, or any other suitable liquid. The quench- 
ing fixture is a \-in. steel plate with 8 holes *%» 
in. in diameter so as to accommodate 8 specimens 
(fig. 9). Metal straps attached to the steel plate 
hold it in position above the quenching bath. After 


' quenching, the specimens are stress relieved at 


250° F to prevent cracking and are then ground 
equally on each, side to 0.394 in. square, thereby 


279 











Fieure 9 Fixture for end-immersion quenching of impact 
specimens 
removing any decarburization that may have 


A hardness survey is then made on one 
face to locate the exact position of the notch. This 
location is marked, a standard V-notch is cut into 
the specimen at this point, and the ends of the 
specimen are cut off at a distance of 1.08 in. on each 
side of the notch. 

During the exploratory work, specimens of a com- 
mercial steel (0.28% C, 1.6% Mn, 0.0015% B) were 
end immersed to different levels, and the hardnesses 
and microstructures at various positions were cor- 
related with those that existed in a standard Jominy 
end-quenched specimen of the same steel. Repeated 
hardness checks and metallographic examinations 
of end-immersed impact specimens showed = con- 
sistently equal hardness and similar microstructures 
over the entire cross section of any selected plane 
parallel to the quenched end. Hardness surveys on 
impact specimens end immersed to various depths 
in brine, together with the Jominy curve for the 
same steel, are shown in figure 10. Examination of 
these curves indicated that any hardness 
ponding to some particular cooling rate) existing in 
the Jominy specimen could be duplicated by end 
immersing the impact specimen. The correlation 
between the normally end-quenched Jominy speci- 
mens and the impact specimens, immersed to diff- 
erent depths in brine, water, and liquid lead-bismuth 
alloy, is shown in figure 11. 

As a result of the exploratory work 
performed on the steel noted, additional tests were 
conducted on one of the experimental steels. An 
impact specimen was end immersed into brine to 
a depth of 7/8 in. and the hardness was measured 
on flats 0.015 and 0.197 in. (approximately center 
below the original surface. The hardness and mi- 
crostructure at both levels checked closely ; the 
curves shown in figure 12 represent the averages of 
the hardness readings plotted against the distances 
from the quenched end. The microstructures of 
both the end-immersed impact specimen and of the 
Jominy specimen of the same steel at 
locations of equal hardness are also shown in this 
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The plotted points represent location f equal hardness and, consequent! 
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bismuth alloy was 300° | 


figure It will be noted that the microstructures at} 
positions of equal hardness are similar and that the 
structure of the impact specimen at any plane of] 
hardness is uniform from just below the surface to 

the center. An additional example of the uniformity In 


of structure that exists from just below the surface} ° pai 
to the center of an end-immersed impact specimen | same | 
of another steel is shown in figure 13 | sired a 
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Uniformity of microstructure from just below the surface to the center and at different distances from the immersed 


end of an impact specimen of steel BR-17 (0.44% C). 


This specimen was end immersed to a depth of 7/8 in, in brine, 


data on the influence of tempering the experimental 
steels, therefore, would be necessary. 

Impact specimens of many of the boron-free steels 
were end immersed in brine to a depth of % in. 
The original surfaces were removed by grinding and 
a hardness determination was made every \¢ in. 
from the quenched end. These specimens were 
then tempered for 1 hr in succession at 300°, 400 
500,° 600,° 700,° 800,° 900,° and 1,000° F and a 
hardness survey was made of the entire length of 
each specimen after each temper. <A similar pro- 
cedure was followed with conventionally quenched 
Jominy bars using a tempering time of 1 hr 
because of the larger cross section. The resulting 
data are presented in figures 14 and 15. Figure 14 
shows the effect of tempering up to 1,000° F on the 
maximum hardness (hardness at \%, in. from the 


quenched end of the Jominy specimen) as measured | 


Etched with 1° nital < SOO. 


on both Jominy and impact specimens. Figure 15 











differs in that it shows the effect of tempering on the | 


planes of originally equal slack-quenched hardnesses 
of pairs of impact specimens of the steels. 

The curves of figure 15 were used as a guide for 
tempering the test specimens of these particular 
steels. These data were sufficient to serve as 4 
general guide, subject to experimental modification, 
for all of the steels. After slack quenching and then 
stress relieving at 250° F, the planes of equal hard- 
ness were located and the specimens were then tem- 
pered for 1 hr at the proper temperature to give the 
desired hardness. The hardness of the same plane 
was checked in each specimen after tempering to 
assure that the proper hardness level had _ been 
reached. Any specimen which failed to develop the 
proper hardness within +1 R, (most were within 
less than +% R,) was not used. 
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4. Impact Tests—Results and Discussion 


4.1. Test Procedure 


All tests were conducted in a Charpy impact test- 
ing machine of 224-ft-lb capacity. The machine con- 
formed to ASTM specifications except for the radius 
of the striking edge which was 0.047 in. instead of 
0.315 in. All specimens conformed to the ASTM 
standard for Charpy V-notch impact specimens 
(0.394 in. square times 2.165 in. long, with a 45° 
V-notch, 0.079 in. deep, and a root radius of 0.010 in). 
The test temperature varied within the range of 

320° F to +300° F. Specimens tested at other 
than room temperature were held at the temperature 
of test in a liquid bath for at least 15 min prior to 
breaking, and the elapsed time between removal 
from the bath and breaking of the specimens was 
about 3 sec. 

Impact tests were made on pairs of steels slack 
quenched at equal cooling rates to identical hardness 
levels. In addition, tests also were made on many 
of these pairs of steels as slack quenched to equal 
hardness levels and subsequently tempered to lower 
but equal hardness levels. Supplementary tests were 
made on some of the same steels as fully hardened 
and then tempered to the same hardness levels as 
the slack-quenched or slack-quenched and tempered 
steels. 

The impact data obtained were plotted on two 
sets of graphs. The first (figs. 16 to 20) compares 
the impact properties of the various pairs of steels 
slack quenched to equal hardness, both with and 
without tempering. Data on the properties of any 
one steel, therefore, often appear in more than one 
graph, depending upon the combination of steels 
involved. In order to assess more directly the effect 
of slack quenching, both with and without tempering, 
upon the impact properties of individual steels, the 
data are replotted in the second set of graphs (figs. 
21 to 27). Where available, data for the fully 
hardened steels are included. 

All fractured specimens were examined to evaluate 
the relative amounts of brittle and ductile fracture. 
It was observed that the gradation in appearance of 
fracture from ductile to brittle was less distinct in 
many slack-quenched specimens than it was in 
fully hardened specimens. This was attributed to 
the mixed structures existing in the slack-quenched 
specimens. Consequently, quantitative determina- 
tions of the relative amounts of brittle and ductile 
fracture in slack-quenched specimens were difficult 
to make in some cases and the values obtained were 
believed to be only approximately correct. For this 
reason the data are not included in the graphs. 


4.2. Slack-Quenched Pairs of Steels 


The test data obtained on pairs of steels, slack 
quenched to the hardness levels indicated, and also 
as slack quenched and tempered, are given in figures 
16 to 20, inclusive. The hardness at which each 
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Impact strengths of slack-quenched pairs of steels 
Code A. 


Figure 16. 


The underlined hardness number in each block indicates the hardness of the 
steels, as quenched. The lower hardnesses were obtained by tempering the 





slack-quenched steels. Maximum hardness: BR-27, 60 R.; BR-22, 58R,.; BR-17, 
55 R.; BR-18, 55 1/2 Ro. 
oo sea c nO oo 4440 408 OF - @044RCc re 
@*- 250 SOK 60 & @* - 2610 30%C +8 ~ erc-e 
Neca Out OC Ht 7“ J - ~~ ae s 2 »c~ 
95% > Yr) : 
om gt -+ on Fn 
> 
d« - 
~ ¥- 
: . ® 
R « 
. 
3 4 
. ’ ~ 
: | 
: 
#0} 
| 
I 
| .f 
Pal 
bao" . 7 2 
_ 
eres ee aw eee aes a ee ee ae ee 





ewrtaatume ° 


Fiat RE 17 Im pact stre naths of slack-quenc hed pairs of steels 
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Code B 
he underlined hardness number in each block indicates the hardness of the 
steels. The lower hardnesses were obtained by tempering the slack-quenched 
steels. Maximum hardness: BR-17, 55 R.; BR-18, 55 1/2 R.; BR-1, 52 R,; 
BR-25, 3 1/2 R.; BR-26, 51 1/2 R.; BR-20, 49 1/2 R 


pair of steels was tested is indicated in the various 
blocks by the underlined value. Each figure gives 
the complete data for a single combination as in- 
dicated by the code thereon and includes, where 
available, a nonboron versus a nonboron steel (left 
column), a nonboron versus a boron steel (center 
column), and a boron versus a boron steel (right 
column). Because of variations in melting, the 
carbon contents often varied slightly from the de- 
sired values. 

A significant observation to be made from these 
figures is that the lower carbon steel of any pair 
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Figure 18 


The underlined hardness number in each block indicates the hardness of the 
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FIGURE 19. 


he underlined hardness number in each block indicates the hardness of the 


steels rhe lower hardnesses were obtained by tempering the slack-quenched 
steels. Maximum hardness: BR-27, 60 R.; BR-28, 60 R.; BR-21, 59 Re; BR-3, 
“1/2 Re; BR-20, 49 1/2 R 


slack quenched to the same hardness invariably 
showed impact properties superior to that of the 
higher carbon steel, regardless of whether or not 
boron was present. With but one or two exceptions, 
where the differences in the impact curves were 
slight, this superiority was maintained after temper- 
ing the slack-quenched pairs. 

The degree of difference in impact properties noted 
between the individual steels of each slack-quenched 
pair was usually greatest in the as-slack-quenched 


284 





assorerc 


enency 


FIGUR! 


The ut 
steels. 1 
steels. \ 
4 Ro; B 


condi 
betwe 
initial 
quence 
carbo! 
greate 
differe 
erally. 
partic 
tures | 
a mor 
fractu 
It n 
of the 
varied 
the hi; 
may | 
shown 
ables 
series 


of the 


43.E 


The 
streng 
shown 
(excep 
split | 
Figure 
heat ( 
steel u 
and as 
son, d 
proper 
tempel 


ad 





ee aauc noe ee -2 46% ” Of - 278 (0 49RC+8 
ee- s% 7 on -26 2 8 ef -30 (0.19% +8 
ace QUENCKHE eon SLACK QUENCHED TO 458 SLACK QUENCHED TO 53M, 
a 
| a3e 
20} 
I 
| ——y 
+ «on 40 & 2% ‘ 
so} 
‘ 
} yr. f 
© ; a A 
<a . 
= so” Lr 
' — - » 
- ae ee a 
i —_— sO  ——— = 
: } 
: son <a (308 sb 30% 
¢ ; > 
. | Sk 
2 
= 4} 
} o—4 
og 
| } 
i p 
fl } *| 
| ‘ 
y oo” 
} ne* “A 
t | oe 
See 
300 -200 OC + +2 ‘ 200 0 + +2 300 -2 10¢ #00 #200 +300 





TEST TEMPERATURE, *F 


FIGURE 20 Impact strengths of slack-quenched pairs of steels 


Code E and Code F 


The underlined hardness number in each block indicates the hardness of the 
steels. The lower hardnesses were obtained by tempering the slack-quenched 
steels. Maximum hardness: BR-28, 60 R.; BR-21, 59 R.; BR-17, 55 Re; BR-30, 
43 R.; BR-24, 4342 Re. 


condition; tempering generally reduced the difference 
between the steels in each pair. The extent of the 
initial difference in impact properties of the as-slack- 
quenched steels was related to the difference in 
carbon content of the steels in each pair; i. e., the 
greater the difference in carbon, the greater the 
difference in impact properties. Tempering, gen- 
erally, raised the impact strength levels of the pairs, 
particularly at the higher and intermediate tempera- 
tures of test. This resulted in the establishment of 
amore pronounced transition from ductile to brittle 
fracture. 

It must be emphasized that, because of the desig1 
of the experimental steels, the total alloy contents 
varied inversely with the carbon contents and that 
the higher amount of alloy in the lower carbon steels 
may be a factor in the better impact properties 
shown by these steels. To separate these two vari- 
ables would have necessitated melting additional 
series of steels not contemplated in the original plan 
of the investigation. 


4.3. Evaluation of Effect of Slack Quenching Upon 
the Impact Properties 


The effect of slack quenching upon the impact 
strength of the steels at certain hardness levels is 
shown in figures 21 to 27, inclusive. Each figure 
(except fig. 27) presents all the impact data on each 
split heat tested, both without and with boron. 
Figure 27 gives data on two steels not from the same 
heat (BR-1 and BR-30). Data are given for each 
steel as slack quenched to the hardness indicated, 
and as slack quenched and tempered. For compari- 
son, data also are given, where available, on the 
properties of the steels as fully hardened and then 
tempered to the required hardness level. 
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The impact properties of the untempered slack- 
quenched steels were lower the higher the hardness 
to which the steels were slack quenched. This was 
attributed to the greater quantity of untempered 
martensite present in the harder steels. Upon tem- 
pering, however, this order was reversed and the 
steels originally slack quenched to the higher hard- 
ness levels had impact properties superior to the 
same steels originally slack quenched to the lower 
hardness levels, thus indicating the deleterious effect 
of slack quenching. This was attributed to the 
larger amounts of tempered martensite present in 
the steels originally slack quenched to the higher 
hardnesses. The steels originally slack quenched to 
the lower hardnesses contained less martensite that 
could be toughened by tempering, and more ferrite 
and pearlite; the two latter constituents are not 
appreciably toughened by tempering. 

As an illustration of the trends just described, 
attention is directed to the data for steel BR-17 
(fig. 23). This steel was slack quenched to three 
levels of hardness, R, 48, 46, and 40. Examination 
of the curves indicates that the impact properties 
improved slightly as the slack-quenched hardness 
decreased. When tempered to 40 and 30 R,, how- 
ever, the impact properties of the steels originally 
slack quenched to 48 R, were better than those slack 
quenched to 46 R,, which, in turn, were better than 
those slack quenched to 40 R,. Another illustration 
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TEST TEMPERATURE, °F 


Ficure 21. Impact strengths of the 0.49%C steels. 

I'he underlined hardness number in each block indicates the hardness of the 
steel. The hardness number on each curve indicates the original slack-quenched 
hardness. Where the underlined hardness is less than the slack-quenched haid- 
ness, the specimens have been tempered. 8S. Q. equals slack quenched. F. } 
equals fully hardened. As fully hardened, both steels had an original hardness 
of 60 Re. 








BR ~ 21(0.46%C, NOB BR-22(0.46%C+8 of this trend is shown in figure 21 for steel BR-28. 




















20 As initially slack quenched to 33 R, and tested at 
40% that hardness, the impact properties were somewhat 
. ra aed superior to those of the same steel initially slack 
oh. quenched to 43 R,, and tested at that hardness 
43 Re * sose However, after tempering both to 30 R,, the steel 
= pe 2.5258 D initially slack quenched to the higher hardness had 
Te aes a [ OO} the better impact properties. When the original 
; oo eh ti slack-quenched hardness was relatively high, thys 
a a indicating only a small amount of nonmartensitie 
© products in the microstructure, the deleterious effeet 
' - . . ¢ ; 
= * 0.408 , of slack quenching upon the impact properties wag 
S x 0—o— yf -* i small. This is illustrated by steel BR—27 (fig. 21) 
5 i y i secsn tested at 40 R.. As slack quenched to 54 R, and 
<2 0 ; then tempered to 40 R,, the steel showed slightly 
: 30 Re —. 28 inferior impact properties only below O°F as com- 
w § . . . e 
Z ; a. pared with the fully hardened (original hardness 
ac 2.48 F equals 60 R,) and tempered steel. The same trend 
was shown by many of the other steels. 
, Faalhiadin Baad Thus, at any particular hardness, the fully har- 
jo dened and tempered steels were usually superior to 
ff the slack quenched and tempered steels. In some 
0 pox o cases, however, they were only equal, or even slightly 
' inferior to the slack-quenched steels over a portion 
" =300 -200 00 +100 +20¢ 300 -200 -100 +100 +2 +3 of the testing range. 
TEST TEMPERATURE, °F The boron and nonboron steels showed one signifi- 
Figure 22. Impact strengths of the 0.46%C steels cant difference in their reaction to tempering in the 
lhe underlined hardness number in each block indicates the hardness of t 
steel. The hardness number on each curve indicates the original slack-quenche 
hardness. Where the underlined hardness is less than the slack-quenched hard 
ness, the specimens have been tempered s. Q equals slack juenched \ test BR-25(0.30%C NOB BR-~26(0.30%C+8 
were made on steels BR-21 and BR-22 as initially fully hardened but int 
condition the hardness was 59 and S58 R., respectively = 
495R _ C sor 7”. ~CO 
4 44% NO 8 BR-i8 44 9 495R 48 
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TEST TEMPERATURE, °F 7 : F 


: : ) 9 > 0", sleels 
FIGURE 23 Impact strengths of the O.4400C steels FIGURE 24 Impact strengths of the 0.30790 els. 


lhe underlined hardness number in each block indicates the hardness of the Che underlined hardness number in each block indicates the hardness of te 
steel. The hardness number on each curve indicates the original slack-quenched steel. The hardness number on each curve indicates the original slack-quenched 
hardness. Where the underlined hardness is less than the slack-quenched hard hardness. Where the underlined hardness is less than the slack-quenched hard: 
ness, the specimens have been tempered. 8S. Q. equals slack quenched. F. H ness, the specimens have been tempered. 8S. Q. equals slack quenched P. H. 
equals fully bardened. As fully hardened, B R-17 had an initial hardness of 55 R equals fully hardened As fully hardened, BR-25 had an original hardness a 
No tests were made on BR-18 as initially fully hardened, but in this condition its 51 R No tests were made on BR-26 as initially fully hardened, but in 8 
hardness was 5542 R condition its hardness was 51! ‘ 
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FiGuRE 25 Impact strengths of the 0 28S %C sleeis 
he underlined hardness number in each block indicates the hardness of the 
stint The hardness number on each curve indicates the original slack-quenched 
hardness. Where the underlined hardness is less than the slaek-quenched hard 
ness, the specimens have been tempered S. Q. equals slack quenched F_H 
equals fully hardened. As fully hardened, both steels had an original hardness of 
92 R 


as-slack-quenched condition. The impact properties 
of the slack-quenched nonboron steels were practi- 
cally always improved by tempering. Whether or 
not the slack-quenched boron steels were improved 
by tempering depended on the differential between 
the hardness as slack quenched and as slack quenched 
and tempered. When this differential was only a few 
points Rockwell C, tempering to 40 R, actually 
caused a slight deterioration in the impact proper- 
ties or, at best, no change. When. however, the dif- 
ferential was large, tempering to either 40 R, or 30 
R,. caused a definite improvement. Examples of this 
phenomenon can be seen in figures 24, 25, 26, and 27. 


5. Summary and Conclusions 


A method was developed and is described for heat- 
treating impact specimens so as to produce uniform 
and accurately controlled slack-quenched microstruc- 
tures and hardnesses at the locations where the 
notches are to be placed. This technique, which 
should prove of considerable use in the study of 
steels, is a modified end-quench procedure. 

A group of steels, based on the 81xx series, varying 
in carbon from 0.19 to 0.49 percent, and with the 
alloy content adjusted at a constant ratio so as to 
raise the hardenability of the lower carbon steels 
and lower the hardenability of the higher carbon 
steels, was melted and fabricated. Each heat was 
split during pouring and boron was added to one- 
half of the melt. Jominy hardenability tests were 
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FIGURE 26 Impact strengths of the 0.19°0C steels 


lhe underlined hardness number in each block indicates the hardness of the 


steel. The hardness number on each curve indicates the original slack-quenched 
hardness, Where the underlined hardness is less than the slack-quenched hard- 
ness, the specimens have been tempered 5s. Y equals slack quenched. F. H. 
equals fullv hardened As fully hardened, BR-23 had an original hardness of 
HR No tests Were made on BR-24 as initially fully hardened, but in this 
condition its hardness was 43! 

made on all steels and the distances from the 


quenched ends at which the hardenability curves of 
various pairs of steels intersected were determined. 
Impact tests were made on Charpy V-notch speci- 
mens of pairs of steels slack quenched to the hardness 
corresponding to the intersection of their respective 
Jominy curves, both as slack quenched and also as 
subsequently tempered. The impact properties of 
some of the slack-quenched steels were compared 
with their fully hardened and tempered counter- 
parts. 

On the basis of the data obtained, the following 
conclusions can be drawn: 

1. When slack quenched at equal cooling rates to 
equal hardnesses and tested at those hardness levels, 
the impact properties of the higher carbon lower 
alloy steel of a pair were always inferior to those of 
the lower carbon higher alloy steel, regardless of 
whether or not boron was present. The greater the 
difference in carbon and alloy content, the greater 
was the difference in impact properties. Tempering 
generally reduced the difference between the impact 
properties and also raised the values. 

2. The hardness level to which any individual 
steel was slack quenched influenced its impact 
properties markedly. As slack quenched only, the 
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FiGuRE 27. Impact strengths of the two steels designated 


The underlined hardness number in each block indicates the hardness of the 
steel. The hardness number on each curve indicates the original slack-quenched 
hardness. Where the underlined hardness is less than the slack-quenched hard 
ness, the specimens have been tempered. 8S. Q. equals slack quenched. F. HH 
equals fully hardened. As fully hardened, BR-1 had an origina! hardness of 
52 Re. No tests were made on BR-30 as initially fully hardened, but in this 
condition its hardness was 43 R 


higher the hardness, the lower were the impact 
properties. Tempering, however, reversed this order 
in that the steels originally slack quenched to a 
higher hardness and then tempered to a lower hard- 
ness had impact properties superior to those of the 
same steels originally slack quenched to a lower 
hardness and then tempered. 

The impact properties of the slack-quenched 
boron steels were improved by tempering to 40 R, 
only when the differential between the hardness as 
slack quenched and as slack quenched and tempered 
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was large. When this differential was small, tem. 
pering actually caused a slight deterioration jp 
impact properties, or, at best, no change. The 
impact properties of all of the slack-quene he d boron | 
a ‘Is were improved upon tempering to 30 R,. 

The impact properties of all slack- -quenched 
ols as tempered to 40 R, were considerably 
improved when the steels were tempered to 30 R,.’ 

At equal hardness levels the impact properties 
of the slack-quenched steels, either tempered or 
untempered, were generally inferior to those of the 
fully hardened and tempered steels. The extent of | 
the deterioration in impact properties caused by 
slack quenching depended on the differential between 
the hardness of the steel as fully hardened and as 
slack quenched. When the tempered slack- -quenched 
steel had an initial hardness close to the maximum | 
possible, the impact properties were but. slightly 
affected. As the initial hardness of the tempered 
slack-quenched steel became lower, the deterioration 
in impact properties became more pronounced. 


The investigation reported in this paper was | 
sponsored by Watertown Arsenal. The authors are 
particularly grateful to Stuart V. Arnold, of Water- 
town Arsenal, whose advice and criticism aided so 
materially in the progress of the work. The authors 
also are deeply appreciative of the conscientious 
assistance of Thomas P. Royston throughout the 
entire course of the investigation and of Nesbit L, 
Carwile, who assisted materially in the latter phases 
of the work. 
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C. The 
a is of C’-Label id’ 
. | Synthesis of C’-Labeled t-Sorbose and 1-Ascorbic Acid 
lenched 
derably 
0R,.° 
thew H. L. Frush and H. S. Isbell 
‘red or 
of the 
‘tent of A method is presented for preparing C-labeled L-ascorbie acid (vitamin C) from 
sed by | labeled barium p-gluconate or p-glucose. The process is adapted from the commercial 
) " synthesis of vitamin C through sorbitol, L-sorbose and potassium di-O-isopropylidene-2- 
etween keto-L-gulonate. In the eyanohydrin synthesis of labeled p-glucose, barium p-gluconate is 
and as used as an intermediate. For preparing labeled L-ascorbie acid, the gluconate is converted 
enched to the lactone. This is efficiently reduced, by means of sodium borohydride, directly to 
xim sorbitol, obviating the intermediate preparation of p-glucose. 
jj um If labeled p-glucose is more readily available than the corresponding gluconic acid (as 
ightly for instance p-glucose-6-C"™ or randomly labeled p-glucose), the sugar is reduced to sorbitol 
npered by the sodium borohydride method; the yield is nearly quantitative. 
ration In a preparation of L-ascorbie acid-6-C™ from barium p-gluconate-1-C™ the yield was 

$1.5 percent, based on the gluconate, or 67.8 percent, based on the intermediate L-sorbose. 

' \ mixture of Methyl Cellosolve, ethylene dichloride, and hydrogen chloride was found to 

| be effective in converting potassium di-O-isopropylidene-2-keto-L-gulonate to L-ascorbic acid. 
r Was | 
TS are 
Vater. | : » . 
led i 1. Introduction | and lactones [7]; for reducing lactones directly to 

‘a SO } . > . — 

athe _ glycitols [8]; and for preparing p-glucose-6-C™ [9]. 
ntioys|  _Carbon-14-labeled L-ascorbiec acid (vitamin C) is | The method reported here for preparing labeled 
it the} useful tool for studying many problems in biology | L-ascorbic acid starts with labeled p-glucose, or 
bit ] and medicine. Before the present work, 1-labeled | barium p-gluconate (which is converted to p-glucono- 
yhases | L-ascorbic acid had been prepared from L-xylosone é-lactone). After reduction of the sugar,’ or lactone, 





lanical 
rought 


through the addition of C'-labeled cyanide [1],? and 
randomly labeled L-ascorbic acid had been prepared 
from randomly labeled p-glucose [2] by the process 
of Reichstein and Griissner [3]. In the present 
commercial process for synthesizing L-ascorbic acid, 
the following substances are successively prepared 





with sodium borohydride, the resulting crude sorbitol 
is not crystallized, but oxidized to L-sorbose with 
Acetobacter suborydans {11, 12]. Preliminary study 
of this oxidation with nonradioactive sorbitol showed 
a maximum yield of L-sorbose of 95 percent in 40 
hr. The yield obtained with crude sorbitol-6-C™, 
however, was lower (69%), and the fermentation 


ion on | from p-glucose: Sorbitol; L-sorbose; di-O-isopropyli- ; 
Trans. dene-L-sorbose ; potassium di-O-isopropvlidene-2-keto- | proceeded _more slowly than that of nonradio- 
Cleve: | L-gulonate; and finally L-ascorbic acid. The process | active sorbitol used as a control. ; 
differs from that of Reichstein and Grissner chiefly | The next intermediate, 2,3-4 ,6-di-O-isopropylidene- 
ity of in the manner of converting L-sorbose to diacetone | L-sorbose-6-C™, was prepared from .-sorbose-6-( , 
tant acetone, and sulfuric acid under conditions previously 





L-sorbose, and potassium di-O-isopropvlidene-2-keto- 
L-gulonate to L-ascorbie acid [4, 5]. Because p-glu- 
cose or D-gluconic acid labeled in position 1, 2, or 6 
is now readily available, it is possible to prepare 
L-ascorbie acid labeled respectively in position 6, 5, 
or 1. This paper reports improved procedures, 
adapted from the commercial process, for preparing 
C-labeled L-ascorbic acid 


2. Discussion of Experimental Methods 


The Bureau has developed efficient methods for 


converting barium carbonate-C™ to sodium cvanide- 
“oS 
C™ [6]; for successively preparing, from sodium 
evanide-C™, C-labeled evanohydrins, aldonic acids, 
See 
Part of & project on the development of methods for the synthesis of radio- 
«tive carbohydrates sponsored by the Division of Research, Atomic Energy 


Commission 


? Figures in brackets indicate the literature references at the end of this paper 


reported by Slobodin [4]. The di-O-isopropylidene- 
L-sorbose, when oxidized with potasstum perman- 
ganate in alkaline solution [3], gave crude potassium 
di-O-isopropylidene-2-keto-L-gulonate in high yield. 
This intermediate, treated with a mixture of hydro- 
chloric acid, methanol, and chloroform, under 
conditions similar to those described by Elger [5], 
vielded about 50 percent of L-ascorbie acid-6-C™. 
However, the reaction mixture turned dark and 
formed a small amount of humuslike material. 
To improve the last step, other solvents and condi- 
tions were studied. It was found that mixtures of 
hydrochloric acid and Methyl Cellosolve (ethylene 
glycol monomethy] ether) considerably improve the 
vield of 1-ascorbie acid, which crystallizes readily 
following the addition of ethylene dichloride. 


—<———___—_—_— 


'Sodium borohydride reduces sugars almost quantitatively to glycitols in 
alkaline solution [10] 
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The accompanying diagram gives a summary of | a test for boric acid with acidified turmeric pape with 
the steps in the synthesis. The over-all vield of | was negative. The material was dissolved in wate} of Fi 
L-ascorbic acid-6-C"™ from barium p-gluconate-1-C™ | and aqueous sodium hydroxide was added until g} itate 
was 41.5 percent, and from t-sorbose-6-C", 67.8 | permanent pink color was obtained with pheng).} The 





percent. phthalein. The solution was then passed through} cont 
Barium p-gluconate-1-C"™ a column containing 30 ml of mixed cation (see fooy} exch 

(89.2) note 5) and anion ® exchange resin. The effluey solu’ 

7, le and wash liquor were concentrated in a vacuum te} abse 
Sorbitol-6-C™ (in solution ) a volume of 20 ml. An assay of the solution pyf cent 
(p-Glucitol-1-( -) a direct count of a portion in formamide [13] showeg} ilu 

| (68.7% >) 1,566 we of activity. Hence, 89.2 percent of the | with 


barium p-gluconate-1-C™ had been reduced, afte} It ¥ 


L-Sorbose-6-C" . P 
- conversion to the lactone. The product, sori} st®? 





84.9! dye . 
(84.9'7/ tol-6-C™, was not crystallized, but used directly fo 18° 
+ j . Sey P + aor ga. (14 , . 
Di-0-isopropvlidene-L-sorbose-6-C"' the preparation of L-sorbose-6-C™. M2 
O4 5°) é 
(94.9'/ 3.2. L-Sorbose-6-C"' mot 
> a : : ; : ‘ re , , . ; L-80 
Potassium di-O-isopropylidene-2- lo study the oxidation of sorbitol to L-sorbose, 1.07 
he to-L-gulonate-6-( solution was prepared from 3.64 g of nonradioactive} .))) 
(84.50) sorbitol, 1.0 g of veast extract, 0.6 g of potassium} — 
. an wes dihvdrogen phosphate and 200 ml of water. . 
L-Ascorbie acid-6-C™ aa. + | | , ; an Bves Ali 
)aeabicgger ipa POPPY quots (25 ml) were placed in 125-ml Erlenmeye 
over-i aT oo /p) ie Te : - : A 
' on ' flasks and sterilized for 15 min at a steam pressure C 
, of 15 lb. Each flask was then inoculated with 5} 4 
3. Experimental Procedures drops of a previously prepared inoculum of Aceto. mill 
bacter suborydans, and incubated at 30° CC. (The} 'U 
3.1. Sorbitol-6-C"' (or D-Glucitol-1-C'') inoculum was prepared by growth, at 30° C for] of f 


15 hr, of a fresh culture of Acetobacter suborydans eal of | 


Sorbitol-6-C"™ was obtained by sodium borohy- | a medium containing 2 g of sorbitol, 0.2 g of yeast} Pe 


dride reduction of p-glucono-d-lactone-1~( a [8]. The | extract, and 40 ml of water.) “rr 
process began with barium p-gluconate-1-C"™ and At intervals, one of the flasks was removed from} 0!" 


was conducted in the following manner: A 1.38-¢ | the incubator, and the volume of solution was} }5 
quantity of barium p-gluconate-1-C'* having an | adjusted to 25 ml at 20° C; 0.2 g of Filtercel (diae} Th 


activity of 1,756 microcuries (uc) * was dissolved in | tomaceous earth) was added, and the mixture was} 





water and passed through a column containing 20 filtered through a drv paper into a 2-dm saeccharim- low 
ml of cation exchange resin.” The combined effluent | eter tube. The optical rotation was measured, and sul 
and washings were concentrated under vacuum in | the yield of sorbose was calculated from the specific | 2 
a 100-ml flask to a thin sirup, which was lactonized | rotation of L-sorbose (—43.4°). Calculated vields} © 
by alternate addition of Methyl Cellosolve and of 23.3. 95.7. 94.4. 92.1. 88.1, and 83.0 percent nes 
evaporation in a current of air. After 8 or 10 | were obtained after incubation times of 20, 40, 45, hot 
evaporations over the course of three weeks, lac- | 49, 66, and 94 hr, respectively. Thus, under the | he] 
tonization was complete and the product appeared | ¢onditions used. L-sorbose reached a maximum in| “=P 
to be entirely crystalline. The flask was then | about 40 hr. fla: 
cooled in an ice bath, and 50 ml of 0.05-\/ aqueous For the preparation of L-sorbose-6-C™, the solu | “* 
boric acid and 10 ml of Amberlite [R 120 resin were | tion described in section 3.1, containing approxi- | &* 
added. With efficient mechanical stirring, 50 ml | mately 4.2 mM of sorbitol-6-C™. and having an Afi 
of freshly prepared 0.3-\7 sodium borohvdride was activity of 1,566 we was divided into two equal parts. | oo 
added during 5 min. Stirring was continued for | Each of these was transferred, with sufficient water | Y® 
30 min, and a second 50-m! portion of the borohy- | to give a volume of 25 ml, to a 125-ml Erlenmeyer |“! 
dride solution was added as before. After an addi- | flask containing 125 mg of yeast extract and 75 mg wh 
tional 30 min, 10 ml of 1-N sodium hydroxide was | of potassium dihydrogen phosphate. The mixture of 
added, and the mixture was stored overnight at | was sterilized at 15 Ib of pressure for 15 min and | di- 
about 5° ©. It was then passed through a column | then incubated at 30° C with five drops of a previ ad 
containing 100 ml of cation exchange resin. The | oysly prepared inoculum of Acetobacter suborydans. - 
effluent and washings were concentrated in a rotating M 


After 72 hr, the fermentation mixtures had approx 
. . . mately the same amount of turbidity as that ob- 

In order to remove boric acid as volatile methyl | served for control experiments with nonradioactive 
borate, 50 ml of methanol was added and the solu- | sorbitol at the point of maximum sorbose content 
tion was evaporated under reduced pressure. This | (49 hr). At. this point the contents of the two 
procedure wis repeated, and after five additions of flasks were combined and treated with 5 ml of de 


vacuum still at 30° C to a volume of about 10 ml. 








methanol, each followed by evaporation to a sirup, 20-percent aqueous zine sulfate, followed by suffi | se 

‘ The work described in sections 3.1-3.2, 3.3, and 3.4-3.5, was carried out wit! cient aqueous barium hvdroxide to give a blue color "? 

materials having specific activities, respectively, of 370, 259, and 60 ye NI SD ¥ lr 
Amberlite IR 120-H, Rohm and Haas Co., Philadelphia, Pa Duolite A4, Chemical Process Co., Redwood City, Calif 
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ctly for 








.-sorbose-6-( it was 


product had an activity of 934 we. An additional 


C pape with bromothymol blue. After the addition of 0.5 g | in 5 ml of water was cooled to 0° C, and a solution of 
1 water} of Filtercel, the solution was filtered, and the precip - | 0.43 g of potassium permanganate in 10 ml of water 
until gf itate was washed thoroughly, and then discarded. | was added dropwise, with shaking, over the course of 
Pheno}.| The filtrate was passed slowly through a column | ¥s hr. Then the mixture was mechanically shaken 
throug,| containing 30 ml of a mixture of cation and anion | for 4 hr. At the end of this time it was heated to 50 
ee foo} exchange resin; a conductivity test of the combined | C to coagulate the manganese dioxide, and 2 ml of 
€ fHlueyt solution and = washings showed nearly complete | ethanol was added to remove unreacted perman- 
uum tf absence of ionic impurities.’ lhe product was con- | ganate. The solution was filtered, and the residue 
tion by centrated under reduced pressure to a thick sirup, | washed with water and discarded; the filtrate was 
showed diluted with about 1 ml of methanol, and treated saturated with carbon dioxide and freeze-dried. To 
of the} with 2-propanol to the point of incipient turbidity. | recover any remaining di-O-isopropylidene-t-sorbose, 
l, afte} It was then seeded with 1-sorbose and allowed to | the dry residue was extracted 4 times with 10-ml 
sorbj.| stand for several days, after which a cropoferystalline | portions of ether. However, the ether extract con- 
separated; the recrystallized | tained less than 1 ye, and hence was discarded. 


The main product, containing 121 yc of carbon-14, 


142 uc of L-sorbose-6-C' was separated from the | was used without purification for the production of 
mother liquor by use of a total of 500 mg of carrier | ascorbic acid, 
} ,-sorbose in three portions. The total radioactivity, 
cell 1,076 we, was 68.7 percent of that of the crude | 3.5, Conversion of Potassium Di-O-isopropylidene-2- 
al te sorbitol-6-C™. keto-L-gulonate-6-C" to L-Ascorbic Acid-6-C™ by 
an Hydrogen Chloride in Methanol 
~ Alief 3.3. Di-O-Isopropylidene-L-sor bose-6-C" —e 
imeye _ ; ah , al ¥ — 
reall One ml of concentrated sulfuric acid was added The solution of crude potassium di-O-isopropyli- 
with ;} dropwise to 25 ml of dry acetone at 0° C. Twenty dene-2-ke to-L-gulonate-6-C™ (2.0 mM_ having an 
Acet.| huilliliters of this solution was placed in a 200-ml activity of 121 we) was passed through a column 
(The round-bottomed flask, together with 540 mg (3mM) | contaiming 10 ml of cation exchange resin. The 
C for of finely powdered L-sorbose-6-C™ having an activity effluent was concentrated under vacuum to a sirup, 
which was transferred with a little ethanol to a 10-ml 


Ans On 





of 777 ue (see footnote 4). The mixture was stop- 
pered, and shaken mechanically until all of the 


glass ampoule. The sirup was then dehydrated by 


F yeast L-sorbose had dissolved (1 to 2 hr). The pale yellow | the repeated addition and evaporation of absolute 
| from} solution was placed first in a refrigerator at 0° C for | ethanol, and finally | was dissolved in 0.2 ml of 
n wast 18 hr, and then in a freezer at —20° C for 24 hr. methanol containing 20 percent of hydrogen chloride. 
| (dig.| The cold material was diluted with a solution of 3 g After the addition of 2 ml of chloroform, the ampoule 
e wast of potassium carbonate in 25 ml of ice water, fol- | Was sealed and stored at room temperature. — In the 
1arim- | lowed by 100 ml of acetone. Precipitated potassium | Course of 3 weeks, a crop of L-ascorbic acid-6-( 
1. and | sulfate was removed by filtration, washed with ace- crystallized from solution; the crystals were sepa- 
pecifie | tone, and discarded. The acetone solution was rated and washed with acetic acid. After recrystal- 
vields | concentrated on a rotary vacuum still nearly to dry- | lization from glacial acetic acid with the addition of 
areent | ness, and the residue was carefully extracted in a chloroform to the point of incipient turbidity, & crop 
0. 43 | hot water bath with successive 50-ml portions of of L-ascorbic acid-6-( having an activity of 71 ye 
1 the | heptane. The extracts were removed by a vacuum | Ws obtained, The dark-colored mother liquors 
were diluted with an equal volume of toluene and.con- 


siphon arrangement to a 200-ml round-bottomed 


“| flask. Following the fourth extraction, the residue | centrated in a rotary vacuum still to a thick sirup. 
soly. | Was dissolved in acetone, and the solvent was By use of 400 mg of L-ascorbic acid as carrier, and 
proxi- | eVaporated as completely as possible under vacuum, | cocrystallization from acetic acid, material having 
ir an | After one additional extraction of the residue, the | 8 additional 20 ue of radioactivity was obtained, 
parts. | combined heptane extracts were concentrated under | © give a radiochemical yield of 75.2 percent. 
water | Yacuum to a sirup, which crystallized when seeded 
never | with di-O-isopropy lidene-L-sorbose. The material, 3.6. Conversion of Potassium Di-O-isopropylidene-2- 
5 mg when recrystallized from hot heptane, had an activity keto-L-gulonate to L-Ascorbic Acid by Hydrogen 
xture | Of 611 we. By use of 500 mg of nonradioactive Chloride in Methyl Cellosolve 
-and | “i-O-isopropylidene-L-sorbose, material having an 
yrevi- additional 49 uc of activity was recovered from the To study the conversion step, l-mM samples of 
dans, | ™other liquor, to give a total radiochemical vield of | potassium di-O-isopropylidene-2-keto-1-gulonate in 
roxi- | *4.9 percent. | Methyl Cellosolve-hydrogen chloride mixtures were 
t ob- ; aoeee stored at 40° C under nitrogen in sealed glass tubes, 
ctive 3.4. Potassium Di-O-isopropylidene-2-keto-L- and the yield of L-ascorbie acid was delrminal by 
itent gulonate-6-C titration with iodine in acid solution. The results 
two A solution containing 500 mg of di-O-isopropyli- | @re reported in table 1. Relatively little color 
il of dene-t-sorbose-6-C™ (having an activity of 128 ye, | developed, and the yields of t-ascorbic acid were 
suffi- | see footnote 4) and 0.225 ¢ of potassium hydroxide | Somewhat higher than those obtained in comparable 
color | ) experiments with the conventional methanol-hydro- 





’ Test made with Barnstead Purity Meter, Barnstead Still and Sterilizer Co., 


Boston, Mass 


gen chloride mixtures. 
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In view of the results given in table 1, the following 
procedure was selected for converting the s/eto-acid 
to t-ascorbie acid: A 2-mM sample of potassium 
di-O-isoprop ylidene-2-keto-L-gulonate in aqueous so- 
lution was passed through a column containing 10 
ml of ice-cold cation exchange resin. The effluent 
was concentrated under vacuum to 10 ml and was 
then freeze-dried. The residue was dissolved in a 
mixture of 1 ml of concentrated hydrochloric acid 
and 7.5 ml of Methyl Cellosolve. The solution was 
diluted with 15 ml of Methyl Cellosolve and stored 
under nitrogen at 40° for 1 week. Then it was con- 
centrated under vacuum to about 1 ml, diluted with 
5 ml of Methyl Cellosolve, reconcentrated to 1 ml, 
and diluted with 5 ml of ethylene dichloride. In the 


course of 24 hr, a crop of crystalline L-ascorbic acid 
The mother liquor was separated, and the 


formed. 
product was washed with acetone. It was recrystal- 
lized by dissolving it in Methyl Cellosolve, concen- 
trating the solution under vacuum to a low volume, 
and adding 5 ml of ethylene dichloride. The re- 


acid f on 


Veth 


TABLE 1 Yields of ascorbic 
propylide ne-2-keto-t-qulonate n 


chloride mixtures 


Rearrangement 1 


Methyl! 
Cellosolve 


* One millimole of potas 
materials listed below 
> Calculated from titrati 


crystallized product weighed 265 mg; titration y 
iodine showed the presence of 1.49 mM of L-aseg 
acid. The mother liquor contained 0.20 mM 
analysis, making the total yield in the last step] 
mM, or 84.5 percent. , 


The authors are grateful to Galen Wampler f 
assistance in the experimental work. 


4. References 


King, Seience 111, 
Burns, and C. G 
1952). 


Gibbs, 


fi] J. J. Burns and C. G. 
L. L. Salomon, J. J 
Chem. Soc. 74, 5161 

[2] a Bothner-Py, M 
Science 112, 363 (1950). 

T. Reichstein and A. Griissner, Helv. 
311 (1934). 

Ya M. Slobodin, J Gen Chem. | 
1947); C. A. 42, 871 (1948). 

I. Elger, Festschr. Emil Barell 1936, 229-37; C. A, 
2585 (1937): U. S. Patent 2,129,317 (1938). 

J.D. Moyer and H.S8. Isbell, Anal. Chem. 29, 393 (195 

H. 8S. Isbell, U. S. Patent 2,606,918 (1952): H. S. Ishe 
J. V. Karabinos, H. L. Frush, N. B. Holt, A. Schwebe 
and T. T. Galkowski, J. Research NBS 48, 163 (195 
RP2301; H. S. Isbell, H. L. Frush, and R. Schaff 
thid 54, 201 1955 RP25S1 

H. L. Frush and H. 8. Isbell, J. Am. Chem. Soe. 
2844 (1956). 

R. Schaffer and H. S 
(1956) RP2667. 

M. Abdel-Akher, J. K. Hamilton, and F. Smith, J. Ag 
Chem. Soc. 73, 4691 (1951); B. Lindberg and | 
Misiorny, Svensk Papperstidn. 55, 13 (1952); C. J 
46, 7942 (1952 P. S. Skell and J. G. Crist, Nat 
173, 401 (1954 

fll] P. A. Wells, J. J. Stubbs, L. B. 
toe, Ind. Eng. Chem. 29, 1385 

112] H. S. Isbell and J. V. Karabinos, J 
138 (1952) RP2334. 

[13] A. Schwebel, H. S. Isbell, and J. D 
NBS 63, 221 (1954) RP2537 


257 (1058 


King, J. J 


and R ( A ndersg ’ 


Chim. Aeta 


S.5.R.), 1%, 


Isbell, J. Research NBS 56, 19 


Lo« kwood, and E. 
(1937 
Research NBS 


Mover, J. Resear 


Wasnineron, April 18, 1957. 

















DIVISION OF PUBLIC DOCUMENTS 
WASHINGTON 25, D. C. 


OFFICIAL BUSINESS 


UNITED STATES PENALTY FOR PRIVATE USE To ap 
‘GOVERNMENT PRINTING OFFICE 


PAYMENT OF POSTAGE, $300 
(ePo0 ; 















- Journal of Research of the National Bureau of Standards Vol. 59, No, 4 





Contents 


Formulas for inverse osculatory interpolation in the complex plane. 


Herbert E. Salzer. RP2792 
An X-ray study of textural stresses in two-phase alloys. C. J. Newton 
and H. C. Vacher. RP2793 
Dissociation of 4-chloro-4’-aminodiphenylsulfone. Elizabeth E. Sager 
and Fleur C. Byers. RP2794 


Speed of sound in water by a direct method. Martin Greenspan and 
Carroll E. Tschiegg. RP2795 


Pattern of a flush-mounted microwave antenna. JamesR. Wait. RP2796 


Standards for pH measurements from 60° to 95° C. Vincent E. Bower 
and Roger G. Bates. RP2797 


A numerical study of Dedekind’s cubic class number formula. Harvey 
Cohn. RP2798 


Impact properties‘of slack-quenched alloy steels. Melvin R. Meyerson 
and Samuel J. Rosenberg. RP2799 


Synthesis of C-labeled L-sorbose and L-ascorbic acid. H. L. Frush 
and H. S. Isbell. RP2800 





[A list of recent NBS publications is given on page 3 of the cover.) 






October 16 





Page 


233 


239 


245 


249 
255 


261 


265 


273 














Subscription price: $4.00 a year; $1.25 additional for foreign mailing. 
UNITED STATES GOVERNMENT PRINTING OFFICE, WASHINGTON : 1957 


Bits bey the Stpertatendeat of Documents, U. S. Government Printing Office, Washington 25, D.C. Price 60 cents (single oc 





